HeanoaMced 
tuA,  tU& 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATIffll 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  p\irpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  -whatsoever;  and  the  fact  that  -the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other- 
■vd.se  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  pennission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


4/^  ASD  TECHNICAL  REPORT  61-192 

lo 


A  THEORETICAL  AND  EXPERIMENTAL  STUDY  OF  GAS  FLOW 

1  THROUGH  CLOTH  OVER  A  RANGE  OF  PRESSURES 

AND  TEMPERATURES 


CO 

CD 

CO 


«£C  cn 

o<c 


FREDERICK  O.  SMETANA 
UNIVERSITY  OF  SOUTHERN  CALIFORNIA 


A 


SEPTEMBER  1961 


•V 


AERONAUTICAL  SYSTEMS  DIVISION 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation,  the 
United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever; 
and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way  supplied 
the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation,  or 
conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Qualified  requesters  may  obtain  copies  of  this  report  from  the  Armed  Services  Tech¬ 
nical  Information  Agency,  (ASTIA),  Arlington  Hall  Station,  Arlington  12,  Virginia. 


This  report  has  been  released  to  the  Office  of  Technical  Services,  U.  S.  Department 
of  Commerce,  Washington  25,  D.  C.,  for  sale  to  the  general  public. 


Copies  of  ASD  Technical  Reports  and  Technical  Notes  should  not  be  returned  to  the 
Aeronautical  Systems  Division  unless  return  is  required  by  security  considerations,  con¬ 
tractual  obligations,  or  notice  on  a  specific  document. 


ASD  TECHNICAL  REPORT  61-192 


A  THEORETICAL  AND  EXPERIMENTAL  STUDY  OF  GAS  FLOW 
THROUGH  CLOTH  OVER  A  RANGE  OF  PRESSURES 

AND  TEMPERATURES 


FREDERICK  O.  SMETANA 

ENGINEERING  CENTER 
UNIVERSITY  OF  SOUTHERN  CALIFORNIA 
LOS  ANGELES  7,  CALIFORNIA 


SEPTEMBER  1961 


DIRECTORATE  OF  MATERIALS  AND  PROCESSES 
CONTRACT  No.  AF  33  (616) -7102 
PROJECT  No.  7320 


AERONAUTICAL  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


900  -  January  1962  -  16-674  S  675 


FOREWORD 


This  report  was  prepared  by  University  of  Southern  California 
under  Contract  No.  AF  33(616)-7102.  This  contract  was  initiated 
under  Project  No.  7320,  "Fibrous  Materials  for  Decelerators  and 
Structures,"  Task  No.  73203,  "Fibrous  Structural  Materials.  "  The 
work  was  administered  under  the  direction  of  Directorate  for  Materials 
and  Processes,  Deputy  for  Technology,  Aeronautical  Systems  Division, 
with  Captain  Clarence  0.  Little,  Jr.  acting  as  project  engineer. 

This  report  covers  work  conducted  from  February  I960  to 
May  1961. 

The  author  wishes  to  thank  Dr.  H.  T.  Yang  for  his  careful  and 
critical  review  of  the  manuscript  and  for  his  assistance  in  clarifying 
the  theoretical  aspects  of  the  flow  phenomena. 

Mr.  Kenneth  Yost  contributed  materially  to  the  experimental 
portion  of  the  program  by  conducting  many  of  the  tests  and  keeping  the 
apparatus  functional. 

The  author  also  wishes  to  thank  the  project  monitor,  Capt. 

C.  O.  Little,  Jr.  for  his  assistance  in  securing  reference  material, 
arranging  conferences  with  users  of  permeability  information,  and 
becoming  thoroughly  acquainted  with  the  problems  of  and  results  of 
the  work.  His  efforts  have  been  responsible  for  widening  the  scope 
of  the  investigation  and,  it  is  hoped,  have  resulted  in  a  more  useful 
work. 


ABSTRACT 


Measurements  were  made  of  the  permeability  of  five  fabrics  at 
downstream  pressures  from  sea  level  to  150,  000  pressure  drops 
across  the  samples  of  1  mm  Hg  to  900  mm  Hg,  and  stagnation  temper¬ 
atures  from  300°K  to  930°K.  In  addition  to  this  basic  information,  the 
investigation  sought  to  provide  a  means  of  predicting  high  altitude  results 
from  those  at  sea  level.  It  was  found  that:  (a)  the  geometry  of  the  test 
apparatus  can  have  a  marked  influence  on  the  results,  (b)  the  major 
elastic  effects  on  permeability  arise  from  the  change  in  fabric  pore 
inclination  with  load  rather  than  through  simple  extension  of  the  yarn 
itself,  (c)  viscous  effects  are  present  for  all  except  the  very  highest 
pressures,  (d)  rarefaction  effects  appear  at  altitudes  above  about 
60,000',  (e)  the  most  satisfactory  model  for  explaining  the  results 
appears  to  be  one  likening  the  flow  to  that  between  two  noninteracting 
cylinders,  and  (f)  as  long  as  the  fabric  retains  its  elasticity  and  does 
not  take  a  permanent  set,  temperature  changes  affect  the  permeability 
only  insofar  as  the  air  density  is  changed.  Some  of  these  findings  have 
not  as  yet  been  described  quantitatively;  for  this  reason,  a  satisfactory 
prediction  procedure  is  not  yet  available. 
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I.  INTRODUCTION 


The  permeability  to  gas  flow  is  one  of  the  physical  properties 
of  cloth  of  most  interest  to  parachute  designers  because  of  the  relation 
which  it  has  to  the  inflation  ability,  stability,  drag,  and  opening  shock 
of  the  parachute.  Permeability  measurements  on  parachute  fabrics 
have  been  made  for  many  years  but  the  realization  that  the  permeability 
will  change  with  pressure  ratio  as  well  as  with  AP  is  of  fairly  recent 
origin.  Since  the  flow  mechanisms  through  cloth  are  still  relatively 
poorly  understood,  despite  intensive  efforts  to  obtain  analytical  repre¬ 
sentations  for  them,  they  must  still  be  investigated  experimentally. 

This  involves  the  use  of  vacuum  facilities  of  considerable  capacity  and 
is  therefore  both  expensive  and  time  consuming.  Hence,  it  was  the 
intent  of  the  present  work  to 

(a)  extend  the  range  of  pressures  over  which  accurate  measure¬ 
ments  of  permeability  on  selected  fabrics  have  been  made,  and 

(b)  study  the  results  of  these  and  other  tests  in  an  effort  to 
describe  at  least  the  laws  by  which  measurements  at  a  single  downstream 
pressure  (atmospheric)  can  be  converted  to  those  at  any  other  downstream 
pressure. 

Results  of  previous  studies  had  indicated  that  the  region  of  chief  interest 
was  that  for  which  AP  <  190  mm  Hg. 


Manuscript  released  May  1961  for  publication  as  an  ASD  Technical 
Report. 
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Since  changes  in  temperature  ciffect  the  density  and  viscosity 
of  the  air  flow  and  may  affect  the  elastic  properties  of  the  cloth,  an 
additional  object  of  the  present  program  was  to  measure  for  the  first 
time  the  effects  of  high  gas  temperature  on  permeability  at  various 
simulated  altitudes.  It  was  desired  also  that  any  effects  of  temperature 
found  in  this  investigation  be  included  in  the  expressions  for  converting 
results  at  one  altitude  to  those  at  another. 

The  materials  investigated  were: 

(a)  Nylon  Fabric,  Mil-C-7350,  Type  I 

(b)  Nylon  Fabric,  Mil-C-8021,  Type  II 

(c)  Glass  Fabric,  Mil-F-9084,  Type  VIII 

(d)  An  experimental  organic  fiber  (HT  -  1)  woven  to  meet 
the  requirements  of  Mil-C-7350,  Type  I 

(e)  A  high  porosity  stainless  steel  fabric.  The  fabric  is 
constructed  from  wire  0.0014”  m  diameter  and  is  a 
broadcloth  with  approximately  220  wires  per  inch  in 
both  warp  and  fill.  The  solidity  is  approximately  52%. 

The  term  permeability  as  used  herein  refers  generally  to  the 
mass  flow  which  will  pass  through  a  unit  area  of  the  material.  This 
differs  somewhat  from  the  more  common  usage  in  which  permeability 
refers  to  the  number  of  standard  cubic  feet  per  minute  which  will  pass 
through  a  square  foot  of  the  material.  The  conversion  from  one  to  the 
other,  however,  involves  only  multiplication  or  division  by  the  actual 
gas  density. 
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II.  REVIEW  OF  PREVIOUS  WORK 


Although  screens  and  cloth  are  similar  in  many  respects  from 
an  aerodynamic  viewpoint,  they  appear  to  have  been  treated  quite 
differently  in  the  literature  on  gas  flow  through  porous  media.  For 
example,  the  screen  data  is  frequently  presented  as  the  variation  of 
AP/q  with  approach  Reynolds  ntimber  or  Mach  number,  while  the  cloth 
data  is  most  often  reported  in  terms  of  the  variation  of  the  volume  rate 
of  flow  through  a  unit  area  with  pressure  drop  (AP).  Possibly  this  dif¬ 
ference  in  presentation  was  a  result  of  the  end  use  of  the  data,  that  on 
screens  having  been  used  primarily  by  wind  tunnel  designers  and  others 
to  calculate  the  power  losses  resulting  from  the  use  of  screens  to 
improve  velocity  profiles.  That  on  cloth  was  used  primarily  by  parachute 
designers  to  correlate  the  performance  of  parachutes.  Nevertheless,  the 
aerodynamic  similarities  present  should  permit  the  use  of  screen  data  in 
the  study  of  the  mechanisms  of  flow  through  cloth.  Screen  data  have  the 
advantage  that  the  geometrical  configuration  involved  is  better  defined 
and  less  subject  to  change  with  air  load  than  is  that  of  cloth.  Conse¬ 
quently,  screens  can  be  analyzed  more  readily  than  cloth.  In  spite  of 
this,  however,  Cornell  (Ref,  1)  concluded  that  losses  through  screens 


were  best  determined  empirically. 


A.  FLOW  THROUGH  PERFORATED  PLATES 


*  In  an  effort  to  reduce  number  of  factors  influencing  the  results 
and  thus  to  render  the  problem  more  tractable  theoretically,  Cornell 
considered  the  flow  through  closely  spaced  holes  in  a  thin  plate.  With 
this  configuration,  the  losses  in  formation  of  the  jets  which  issue  from 
the  holes  and  those  in  the  plate  are  negligible  if  the  flow  velocities  and/or 
densities  are  high  enough  so  that  the  Reynolds  number  based  on  local 
conditions  in  the  holes  exceeds  10^.  The  entire  loss  can  therefore  be 


assumed  to  take  place  downstream  of  the  plate  as  a  result  of  mixing  of 
the  jets  and  wakes.  This  is  a  problem  which  has  been  studied  to  some 
extent  and  fairly  reliable  theoretical  solutions  are  known.  Following 
earlier  treatments,  Cornell  expressed  the  pressure  drop  across  the 


plate  in  incompressible  flow  as 

Ap  ^  f 

Ir  <^(I-S) 


(1) 


where  <f>  is  the  contraction  ratio  of  the  jet  leaving  the  orifice  and  s 
is  the  ratio  of  the  blocked  area  to  the  total  area  of  the  plate.  They  used 
the  value  of  0(s)  calculated  by  vonMises  in  1917  which  was  confirmed 
by  later  experiment  on  large,  single  orifices.  For  plates  with  large  holes 
and/or  low  solidity  this  equation  appears  to  give  satisfactory  results.  It 
may  be  converted  to  a  form  more  familiar  to  parachute  designers  as 
follows: 
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HL 

A 


r  <^(i-s)  1  /2Ap(P2+ApT 
Ll-(^(l-S)J  ^  ZTx 


(2) 


Some  study  of  the  problem  will  reveal  that  no  matter  what  the 
loss  mechanism  downstream  of  the  plate  —  whether  is  fully  recov¬ 
ered,  whether  all  the  dynamic  head  is  lost,  or  whether  some  condition 
in  between  actually  occurs  —  the  form  of  Eq.  (2)  remains  the  same.  Only 
the  term  in  brackets  will  change.  As  long  as  the  flow  is  incompressible 
it  depends  only  on  the  plate  geometry.  To  illustrate  this  point,  it  is 
rather  easy  to  show  that  if  one  assumes  that  all  the  dynamic  head  is 
lost,  i.  e. ,  that  the  pressure  far  downstream  of  the  plate  is  the  same  as 
the  minimum  pressure  in  the  jets,  then  Eq.  (2)  would  be  written 

m  <^(I-S)  /~2AP(P2+AP) 


So  that  the  conditions  under  which  Eq.  (1)  was  derived  are  not 
violated,  one  must  restrict  its  application  to  cases  where  P2  is  large 
and  AP  small.  For  these  conditions,  Eq.  (2)  could  be  written 


^ _  r  <^(i-s)  -j  /  gPzAp 
A  “  L  |-<^(|-S) 


(3) 


Cornell  attempted  to  extend  the  range  of  applicability  of  Eq.  (1) 
by  including  compressibility  effects  on  the  left  hand  side  of  the  equation 
and  using  empirically  determined  values  for  <^(s).  These  were  for 
^  2./^  I  ^  0.528,  however.  The  agreement  which  Morgan  (Ref.  2) 
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obtained  between  this  revised  equation  and  his  experimental  results 
was  increasingly  poor  for  ^2/^1  ^  Cornell's  relation  seemed 

to  consistently  underestimate  the  pressure  drop,  indicating  possibly 
that  more  of  the  dynamic  head  was  being  lost  (probably  through  shock 
waves)  than  for  P2/P1  — *  1*0  that  <^(s)  at  the  pressure  ratios 
involved  was  different  than  assumed.  Despite  this  lack  of  agreement 
it  seems  clear  from  fundamental  considerations  that  one  can  always 
represent  the  mass  flow  through  porous  plates  at  high  Reynolds  numbers 
by  the  functional  form 

-5-=  f  (s,^p,)g(P2/P,)  (4) 

where  g(^2/^0  familiar  isentropic  nozzle  flow  function.  Whether 

one  can  arrive  at  a  satisfactory  prediction  for  f(s,  P^/P^),  however,  is 
another  matter.  In  order  to  do  this  conveniently,  one  must  be  able  to 
develop  an  expression  for  the  mixing  losses  which  is  simultaneously 
valid  for  all  of  the  following  cases; 

(1)  Wherein  the  individual  jets  are  perturbed  primarily  by  wakes 
such  as  occurs  when  the  solidity  is  less  than  about  0.4  (Ref.  3). 

(2)  Wherein  the  individual  jets  —  because  of  relatively  large 
distance  between  them,  as  for  example,  in  a  plate  with  high  solidity  — 
react  rather  violently  upon  contact. 

(3)  A  situation  between  cases  (1)  and  (2). 


(4)  Wherein  the  pressure  ratio  across  the  plate  is  high  enough 
to  produce  expansion  and  acceleration  of  the  jets  (see  Appendix  I  for 
an  elaboration  of  this  point). 

(5)  Wherein  the  jet  flow  may  be  assumed  to  be  incompressible 
This  is  obviously  a  difficult  task,  one  which  has  as  yet  not  been  com¬ 
pleted  successfully.  One  could,  perhaps,  attempt  to  split  the  task  into 
parts,  i.e.  ,  deVelop  expressions  valid  only  for  given  regions  of  ^2/^! 
and  s,  but  this  then  leads  to  problems  in  determining  the  applicability 
of  the  results  and  fairing  between  regions  of  applicability.  At  the 
present,  therefore,  one  must  conclude  that  while  exact  values  are 
available  for  certain  special  cases  of  f(s,  P2/P2)  and  the  qualitative 
behavior  is  understood  in  general,  no  general  quantitative  description 
is  within  sight.  Thus,  one  could  not  be  expected  to  predict  accurately 
the  mass  flow  through  an  arbitrary  perforated  plate  for  any  arbitrary 

AP  and  P2  even  if  the  results  are  known  to  be  insensitive  to  changes 
in  Reynolds  number.  Similarly,  it  is  virtually  impossible  to  predict 
the  geometric  solidity  of  a  plate  from  experimental  data  except  for 
isolated  cases. 

B.  FLOW  THROUGH  ROUND  WIRE  SCREENS 

In  reviewing  earlier  work  on  incompressible  flow  through  screens, 
Cornell  found  that  the  semi-empirical  equation  proposed  by  Wieghardt 
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0 


AP 


6S 

Re''^{l-S)2 


for  60  <  Rg  <  1000 


and  that  proposed  by  MacDougall 


Ap 


-1.27 

33,93  S(l-S) 

Re  I  +  ( I  -  S )  -5 


for  .006  <  Rg  <  20, 


•where 


Re 


U.d' 

2/(|-S) 


(5) 


(6) 


(7) 


and  d'  is  the  ■wire  diameter,  agreed  well  with  his  own  results  and  that 
of  others.  For  1000  <  Rg  <  17,000,  Cornell  found  that  the  results 

were  for  practical  purposes  independent  of  Reynolds  number. 

Apparently  the  initial  form  of  these  results  were  obtained  by 
considering  the  screen  to  be  made  up  of  a  number  of  independent 
circular  cylinders  and  converting  empirical  Cq  ~  Re  date  to  a 
suitable  form.  Since 


D  Ri  APA, 


Cq  RJ 


Ap 


(8) 

(9) 


On  this  basis,  one  would  assume  that 


Ap 


Re^O.75 


for  Rg  <  2 


(10) 


Ap 


Red.  0.25 


for  10  <  Rg^  <  1000  (11) 
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q 


constant 


for  10^  <  Re  ,  <  2  x  10  (12) 

d 

Little  or  no  data  is  available  to  indicate  the  behavior  of  C£)  when 
M  >  0.4  and  Re^j  ^  10^.  One  would  speculate,  however,  that  a 

Reynolds  number  effect,  superimposed  on  the  usual  Mach  number  effects 
is  to  be  expected.  A  more  quantitative  analysis  of  screen  flow  based  on 
this  approach  is  given  in  Appendix  I. 

Study  of  Eqs.  (5)  and  (6)  reveals  that  they  can  be  written  in  the 

form 

=  h(s,  P2/P1.  Re)  g(R2/Pl)  (^3) 

The  losses  which  occur  upstream  of  and  in  the  screen  can  be  included 
in  the  dependence  on  Reynolds  number.  (This  question  is  discussed  at 
some  length  in  Appendix  I.  )  Thus,  the  form  is  basically  the  same  as 
for  Eq.  (4)  except  for  the  inclusion  of  the  Reynolds  number  dependence. 


C.  FLOW  THROUGH  CLOTH 

With  this  background,  one  can  proceed  somewhat  more  know- 
ledgably  to  review  the  efforts  made  to  predict  the  flow  through  cloth. 
During  the  discussion  which  follows  it  should  be  remembered  that 


compared  with  screens  or  perforated  plates,  most  cloth  samples  are 


quite  elastic  and  have  a  relatively  high  solidity.  The  high  solidity 
usually  leads  to  complete  loss  of  the  dynamic  head  -while  the  elasticity 
leads  to  changes  in  effective  solidity  with  AP. 

Lavier  (Ref.  4)  after  studying  earlier  work  on  the  flow  of  air 
through  cloth  and  results  of  his  own  experiments  concluded  that; 

(a)  For  very  low  pressure  differentials  (i.  e.  ,  <50  iTi.H20) 

across  the  cloth,  the  flow  "is  like  that  of  viscous  flow  through  a  porous 
medium”.  j^Study  of  works  on  this  subject,  for  example  Ref.  5,  shows 
that  for  such  cases  rii/A  ~  A  P.  A  similar  result  is  observed  in 

Eq.  (6)]  . 

(b)  For  higher  pressure  differentials,  "the  air  flow  through  the 
cloth  may  be  like  the  air  flow  through  a  screen  or  grid”.  Presumably, 
the  dependence  of  the  mass  flow  through  the  cloth  on  Reynolds  number 
is  different  than  in  (a). 

(c)  For  still  higher  pressure  differences,  "compressibility 
effects  will  be  experienced”,  as  well  as  elastic  effects. 

(d)  Specification  of  the  proper  characteristic  length  presents 
the  chief  problem  in  obtaining  a  better  description  of  the  viscous  effects. 

Lavier  also  reviewed  efforts  to  apply  the  theory  of  Green  and 
Duez  (Ref.  6)  to  the  flow  through  cloth.  In  this  theory,  it  is  assumed 
that  the  flow  is  incompressible  and  isothermal  and  that  the  pressure 
drop  across  the  cloth  can  be  represented  as  a  linear  sum  of  a  viscous 
term  (for  which  the  pressure  drop  is  linearly  proportional  to  flow 
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velocity)  and  an  inviscid  term  (for  which'the  pressure  drop  is  propor¬ 
tional  to  the  square  of  the  velocity);  the  coefficients  are  left  undetermined 
to  be  evaluated  from  experimental  data  and  are  presumably  functions  of 
the  geometry  only.  He  found  that,  generally,  for  Ap  >  50  in.H20»  the 

theory  predicted  increasingly  less  mass  flow  (compared  with  the 
measured  values)  with  increasing  pressure  drops  for  both  parachute  cloth 
and  100  mesh  screen;  he  was,  therefore,  unable  to  determine  constant 
values  for  the  coefficients.  He  attributed  this  to  disappearance  of  the 
viscous  effects  more  rapidly  than  predicted  by  the  theory. 

In  a  later  report,  Lavier  and  Boteler  (Ref.  7)  concluded  that  the 
flow  through  cloth  for  large  pressure  ratios  across  the  cloth  should  be 
treated  by  a  compressible  theory;  however,  in  calculating  an  effective 
discharge  coefficient,  i.e.,  the  ratio  of  the  actual  mass  flow  to  that 
predicted  by  theory,  they  employed  incompressible  theory  and  based 
their  Reynolds  number  on  the  conditions  in  front  of  the  fabric. 

The  MIT  group  which  examined  these  results,  concluded  (Ref .  8) 
that  such  assumptions  are  tantamount  to  assuming  that  the  flow  through 
the  cloth  is  similar  to  the  flow  through  an  orifice  (presumably  in  incom¬ 
pressible  flow).  This,  they  suggest,  is  incorrect;  rather,  the  flow  can 
be  represented  more  satisfactorily  and  more  correctly  by  assuming  it 
analogous  to  flow  in  an  isentropic  nozzle.  They  attribute  the  pressure 
drop  experienced  across  the  cloth  entirely  to  mixing  losses  downstream 
of  the  fabric.  Implicit  in  their  approach,  although  not  expressed,  is  the 


assumption  that  all  the  velocity  head  is  lost  in  the  mixing  process.  They 
apparently  felt  that  although  the  specification  of  effective  open  area  was 
still  a  procedure  of  uncertain  results,  a  method  whereby  measured 
permeabilities  could  be  extrapolated  to  other  values  of  P2  would  be  very 
useful  in  reducing  the  number  of  tests  required  to  describe  the  perme¬ 
ability  completely.  By  this  device,  the  elastic  effects  need  not  be  known 
explicitly,  and  the  problem  is  reduced  to  describing  the  aerodynamic 
behavior  of  the  cloth  for  a  constant  geometry. 

The  scaling  relations  which  they  developed  are  based  on  the 
variation  in  mass  flow  through  an  isentropic  nozzle.  If  the  elastic 
properties  are  functions  of  AP  alone,  then  according  to  this  analogy, 
the  permeability  at  one  P2  and  AP  may  be  transformed  to  that  at 
another  P2  by  multiplying  by  the  appropriate  ratio  of  g2(P2/Pi)/ 
gl{I>2_/'P])  for  the  same  AP.  They  suggest  that  isentropic  theory  is 
valid  for  AP  >  100  in,H20.. 

The  agreement  with  experiment  achieved  by  these  scaling 
relations  is  only  fair.  In  particular,  the  slopes  of  the  individual 
curves  appear  to  be  different  from  those  given  by  theory  although  the 
magnitudes  of  the  higher  altitude  data  have  about  the  predicted  relation 
to  the  sea  level  results  in  most  cases.  It  would  seem,  therefore,  that 
while  this  approach  may  offer  a  coarse  means  of  extrapolating  sea  level 
data  to  higher  altitudes,  it  is  not  sufficiently  inclusive  or  it  does  not 
employ  a  sufficiently  correct  model  to  permit  accurate  calculation. 
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D.  SPECIFICATION  OF  THE  EFFECTIVE  OPEN  AREA 


The  most  complete  study  to  date  of  means  for  specifying  the 
effective  open  area  seems  to  have  been  given  in  Ref.  9-  In  this  work 
they  define  the  discharge  coefficient  by 


where 


_ m 


2P,  AP 
7^T| 


Af  -  the  open  area  when  AP  =  0  as  determined  from  light 
penetrability  measurements 

Aq  the  ratio  of  the  actual  open  area  which  is  inclined  to  the 
Af  ~  plane  of  the  cloth  to  the  projected  open  area,  i.  e. ,  the 
open  area  as  seen  from  light  penetrability  tests 


(i 


— )  r_L_ 

^  w  '  L  2  cos  0 


1 

^  4  tan  0 


/  1  +  sin  0 
M  -  sin  0 


r  =  the  thread  radius 

w  =  the  projected  spacing  between  threads 

0  =  the  inclination  relative  to  the  plane  of  the  cloth  of  the 
line  between  thread  centers,  i.  e.  ,  pore  inclination 

■^f' 

-  =  the  open  area  increase  due  to  air  load 

Af 


The  pore  Reynolds  number  is  defined  by 


V  r  w(l4  sec9)+  2r(sec0  -  1)~[ 


0  is  determined  from  the  empirical  relation 


6  =  tan“^ 


0.136  v^percent  crimp  in  filling  yarn 


warp  threads  per  inch  x  spacing  of  fill  threads  ^^/thread) 


(17 


r  is  determined  by  measurement  or  from  experimental  relations 
between  yarn  diameter  and  twist,  w  is  then  determined  from  the 
relation  between  yarn  diameter  and  threads  per  inch.  The  change  in 
open  area  due  to  air  load  was  found  from  measurement  of  linear 
extension  and  light  penetrability.  These  data  indicate  that 


-  «  1  +  0.2 


(18) 


for  the  fabrics  tested. 

Assuming  these  relationships  to  be  correct,  it  was  found  that 

Cde  =  Cdg  ^  0  -0.0041  0  (19) 

(0  in  degrees) 

Cd0  =  0  =  (20) 

for  100  <  Re  <  300 

Good  agreement  with  experiment  was  obtained  for  a  number  of  different 
fabrics  for  AP  <  30  "  H^O  and  P2  =  406"  H^O.  Similar  correlation 
for  data  with  AP  >  30  "  H^O  was  not  attempted. 

Several  interesting  conclusions  can  be  drawn  from  this  work: 

1.  From  the  measurements  of  linear  extension  of  the  cloth 
fibers  under  load  it  seems  apparent  that  the  change  in  total 


surface  area  (and  thus,  the  change  in  pore  size  related  to  the 


change  in  area)  is  small,  being  between  7  and  11  percent  f<^r 
the  fabrics  of  interest  in  the  current  study  at  AP  =  680  mm. 

2.  The  change  in  light  penetrability  with  load,  on  the  other 
hand,  is  much  larger  [see  Eq.  (18),  for  example]  indicating 
that  changes  in  6  resulting  from  air  load  are  the  principal 
reasons  for  changes  in  the  effective  open  area.  This  would 
appear  to  be  an  area  for  fruitful  analytical  work  since  the 
momentum  transport  and  thus  the  load  (AP)  on  the  cloth  for 

a  given  mass  flow  depend  on  the  changes  in  direction  which  the 
flow  must  make  in  traversing  the  cloth.  It  is  interesting  also, 
that  the  light  penetrability  measurements  reported  by  Klein, 
et  al,  appear  to  correlate  well  with  the  measurements  of  6  as 
a  function  of  AP  reported  in  Ref.  8. 

3.  The  Reynolds  number  effect  shown  by  Eq.  (20)  seems  to 
be  confined  to  too  small  a  range  to  permit  attaching  much 
credence  to  it.  Generally,  one  would  expect  that  the  Reynolds 
number  dependance  would  be  at  least  as  strong  as  that  for 
turbulent  pipe  flow,  i.e.,  Rg^-^^  to  R^  as  Rg  increases. 

In  the  next  section,  an  attempt  is  made  to  combine  the  signifi¬ 
cant  features  of  the  various  analyses  and  experimental  results  cited 
above  to  obtain  a  more  general  semi- empirical  description  of  gas  flow 
through  cloth. 
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III.  GENERALIZED  ANALYSIS  OF  FLOW 
THROUGH  CLOTH 

It  is  evident  from  the  foregoing  that  a  prediction  of  the  flow 
through  a  given  sample  must  be  based  on  rather  complete  measurements 
of  the  construction  of  the  cloth  and  its  variation  with  load.  It  is  not 
the  purpose  of  this  section  to  deal  explicitly  with  this  portion  of  the 
problem,  however,  and  thus,  no  general  prediction  procedure  can 
results  What  is  sought  is  a  method  whereby  experimental  results  on  a 
group  of  fabrics  can  be  correlated  to  determine  the  influence  of  air 
density.  If  this  is  successful,  then  the  technique  can  be  used  with  fair 
assurance  of  success  to  predict  the  change  in  permeability  with  altitude 
for  a  variety  of  fabrics  from  their  sea  level  permeability  measurements. 
To  do  this,  it  is  necessary  to  know  the  effect  of  elasticity  on  the  perme¬ 
ability  or  to  determine  the  viscous  effects  independently.  Subsection  A 
is  a  study  of  the  latter.  Since  success  in  this  endeavor  would  obviate 
the  need  for  detailed  study  of  the  elastic  effects,  the  analysis  of  these 
effects  —  described  in  Subsections  B  and  C  —  sought  qualitative  results 
only.  It  should  be  pointed  out  in  this  connection  that  even  a  completely 
successful  description  of  the  effects  of  elasticity  on  permeability  by  itself 
would  not  permit  attainment  of  the  objectives  of  the  program,  except  in  a 
wholly  empirical  fashion.  A  satisfactory  physical  description  of  the  flow 
phenomena  is  still  required  to  obtain  an  analytical  expression  valid  over 
the  entire  range  of  interest. 
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A.  EFFECTS  OF  VISCOSITY 


Following  the  approach  used  by  the  MIT  group,  initial  efforts  of 
the  present  study  centered  on  obtaining  a  satisfactory  description  of  the 
mass  transfer  characteristics  of  a  nozzle  when  viscous  effects  are 
important.  It  was  soon  found  that  nothing  short  of  a  numerical  solution 
of  the  general  equations  of  motion  would  yield  analytical  results  in  which 
much  confidence  could  be  placed.  Since  the  proper  boundary  conditions 
for  the  flow  through  cloth  are  at  best  poorly  known,  the  correspondence 
of  a  solution  of  the  nozzle  flow  equations  to  the  flow  through  cloth  may 
not  be  good.  The  very  large  effort  required  to  obtain  numerical  solutions 
was  therefore  not  justified.  As  a  result,  the  effort  was  made  to  describe 
nozzle  flow  in  a  relatively  simple,  semi- empirical  fashion.  The  details 
of  this  analysis  are  given  in  Appendix  A  and  the  results  shown  in  Figs.  Al, 
A2  and  A3.  Basically,  this  procedure  uses  one -dimensional  flow  theory 
modified  by  a  discharge  coefficient.  It  was  found  that  as  long  as  the 
Knudsen  number  of  the  throat  is  less  than  0.01,  and  the  entrance  radius 
is  greater  than  0.25  times  the  throat  diameter,  the  discharge  coefficient 
can  be  expressed  as  a  function  of  Reynolds  number  which  contains  only 
one  empirical  constant;  further,  fairly  crude  theoretical  estimates  of  this 
constant  for  a  particular  geometry  indicate  good  agreement  with  exper¬ 
iment,  particularly  for  Re^j  >  500.  By  writing  the  Reynolds  number  in 

the  form 
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(21) 


I 

I 

i 


iT/4d|i 

where  d  is  the  nozzle  throat  diameter  one  may  then  find  an  explicit 

iTi 

expression  for  the  mass  flow  in  terms  of  Pj,  P21  Tg,  A^i  A]^/A2>  M- 1 
and  the  empirical  constant  B.  This  equation  is  graphed  in  Figs.  1(a) 
and  1(b)  for  two  values  of  d  and  B  =  4.0.  The  effects  of  viscosity  are 
evident  from  a  comparison  of  the  curves  for  the  two  diameters. 

To  obtain  some  feeling  for  the  change  in  appearance  of  the  m  vs. 
AP  curves  resulting  from  a  change  in  the  flow  model,  similar  curves 
were  computed  for  a  sharp-edged  orifice  and  for  an  idealized  screen. 
These  are  shown  in  Figs.  2(a),  2(b)  and  3.  The  curves  for  an  orifice 
were  computed  assuming  that  the  viscous  and  inviscid  effects  could  be 
treated  separately  as  indicated  by  the  data  discussed  in  Appendix  I.  The 
curves  for  the  screen  were  computed  from  published  values  of  Cj;)  vs.  M 
and  Cj)  vs.  for  a  cylinder  according  to  Eq.  (A26).  It  is  evident  upon 
examination  that  there  is  considerable  difference  in  the  magnitudes  and 
slopes  of  the  curves  for  the  three  models,  particularly  at  the  lowe  r 
densities  and  flow  velocities.  Thus,  unless  some  single  model  can  be 
shown  to  represent  flow  through  cloth  it  is  not  possible  to  determine  the 
viscous  effects  independently  of  the  elastic  effects.  Unfortunately,  data 
available  in  the  literature  at  this  time  is  either  too  limited  in  range, 
contains  effects  of  the  test  setup  geometry  (see  Section  V),  or  does  not 
provide  corresponding  results  on  Af/Af  vs,  AP  to  allow  one  to  assign 
a  definite  model. 
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It  will  be  noted,  however,  that  for  small  ranges  of  AP  and 


Ts  all  the  curves  can  be  represented  by  the  relation 


where  G  is  a  constant  which  depends  on  geometry  and  P^,  and  x  is 
a  constant  whose  value  lies  between  0.5  and  1.5.  *  For  inviscid  flow 
X  is  0.5  when  AP  is  small.  When  P^/Pi  <  0.528,  X  increases 

toward  1.0.  For  fully  developed  viscous  flow,  0.5  is  added  to  the  inviscid 
result.  As  the  Reynolds  number  increases,  something  less  than  0,5  is 
added  to  the  inviscid  result. 

Eq.  (22)  being  consistent  with  Eq.  (13),  has  therefore,  quite  general 
validity  if  G  is  also  allowed  to  vary  with  changes  in  loss  mechanism.  It 
will  be  seen  that  Eq.  (22)  can  be  used  to  represent  flow  through  cloth  as 
well  if  X  is  also  made  to  include  elastic  effects.  The  change  which 
temperature  variations  produce  in  the  viscous  and  elastic  characteristics 
can  also  be  accommodated  by  changes  in  x.  Thus,  it  would  appear  that 
the  permeability  characteristics  of  a  particular  cloth  sample  can  be 
represented  quite  adequately  by  evaluating  G  as  a  function  of  pressure 
ratio  and  x  as  a  function  of  m/A,  P2/P1,  and  Tg.  While  these  functions 


The  observation  that  Eq.  (22)  can  be  fit,  using  constant  values  of  x,  to 
all  low  AP  permeability  data  seems  to  have  been  made  first  in  Ref.  10. 
Apparently,  subsequent  investigators  make  no  mention  of  it. 


5,  must,  for  the  present,  be  empirical  they  provide  a  common  basis  for 

a  comparison  of  various  fabrics  which  is  more  indicative  of  the  significant 

aerodynamic  parameters  than  representations  used  heretofore. 

Despite  the  handicaps  imposed  by  such  a  semi-empirical  analysis, 
it  is,  nevertheless,  possible  to  arrive  at  some  conclusions  regarding 
the  significance  of  experimental  data: 

(1)  X  can  be  less  than  0.5  only  if  the  discharge  coefficient 
decreases  with  increasing  Reynolds  number.  This  can  occur  over 
small  ranges  of  Rg  for  which  there  is  flow  separation  from  the  surface 
such  as  through  orifice  plates.  This  would  not  be  expected  for  a  round 
wire  screen  of  fairly  close  mesh  and  for  cloth;  therefore,  any  experi¬ 
mental  data  evidencing  such  a  trend  should  be  regarded  with  extreme 
suspicion. 

(2)  If  X  is  greater  than  1.5,  significant  elastic  effects  are  present. 

(3)  If  both  viscous  and  elastic  effects  are  present,  it  is  extremely 
difficult  to  deduce  their  individual  contributions  to  x. 

Some  elaboration  of  this  last  point  is  in  order.  Presumably,  the 
elastic  properties  at  a  given  AP  and  Tg  are  independent  of  P2;  thus,  it 
should  be  possible  to  find  the  variation  in  x  with  ih/A  at  constant  AP 
if  some  assumptions  regarding  G  can  be  made.  It  is  necessary  to 
assume  that  the  effective  open  area  remains  constant,  and  that  the  vari¬ 
ation  in  G  with  P2  is  directly  related  to  that  for  an  isentropic  nozzle. 
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i.e.  ,  there  is  no  change  in  loss  mechanism  with  P2.  This  procedure 
permits  one  to  determine  the  change  in  the  viscous  contribution  to  x 
with  P2,  but  does  not  allow  one  to  find  the  absolute  value  of  this  contri¬ 
bution  to  X  at  any  P2.  At  some  larger  AP,  the  effective  fabric  pore 
size  increases.  In  order  to  determine  the  change  in  x  due  to  elastic 
effects,  it  is  necessary  to  hold  the  Reynolds  number  constant.  To  do 
this,  however,  it  is  necessary  to  hold  the  ratio  of  m/A  to  the  effective 
pore  diameter  constant.  This  is  obviously  not  possible  so  that  the  best 
one  can  do  is  to  maintain  m/A  constant  as  AP  is  increased  and  hope 
that  the  change  in  dg^^  is  not  sufficient  to  produce  marked  changes  in 
the  viscous  effects  or  loss  mechanism.  If  this  is  true,  then  the  elastic 
contribution  to  x  should  be  the  same  at  all  m/A.  By  iterating  back 
and  forth  until  a  consistent  set  of  values  is  obtained  one  can  probably 
arrive  at  a  reasonably  good  picture  of  the  combined  effects  of  fabric 
elasticity  and  flow  losses,  provided  the  assumptions  in  regard  to  loss 
mechanisms  are  justified.  It  is  evident,  however,  that  this  entire 
process  is  time  consuming  and  because  of  data  scatter  and  lack  of 
knowledge  of  constancy  of  the  loss  mechanisms,  may  be  of  limited 
accuracy.  The  results  should,  nevertheless,  be  valuable  in  making 
qualitative  comparisons  between  various  fabrics  and  establishing 
trends  with  altitude,  load  and  temperature. 
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In  an  effort  to  remove  some  of  the  uncertainty  from  the  afore¬ 
mentioned  procedure,  an  attempt  was  made  to  find  an  analytical,  if 
somewhat  idealized,  description  for  the  change  in  pore  dimension  with 
load  and  its  effect  on  mass  flow.  While  a  complete  description  of  the 
change  in  pore  size  with  load,  rate  of  load,  load  history,  etc.  ,  would 
be  extremely  complex  and  difficult  to  use  even  if  it  could  be  found,  an 
analysis  assuming  the  material  to  be  perfectly  elastic  probably  accounts 
for  most  of  the  change.  Such  an  analysis  can  be  divided  into  two  parts 
for  simplicity. 

B.  CHANGE  IN  SAMPLE  OPEN  AREA  DUE  TO  SAMPLE  DEFLECTION 
WITH  LOAD 

If  one  assumes  that  the  test  sample  deflects  under  load  to  form 
a  portion  of  a  sphere,  then  the  total  surface  area  is 

7 

A  =  it{h^  +  — )  "  (23) 

where  h  is  the  deflection.  If  the  area  closed  to  air  flow  remains 
essentially  constant,  tor  “ecsonable  values  of  h  and  the  deflection  due 
to  load  is  analogous  to  that  of  a  thin  membrane,  then 


and 


"  Z  =  1  +  C  Ap2/3  (26) 

Z  is  the  ratio  of  strained  to  unstrained  fiber  length.  Examination  of 
the  data  in  Ref.  8  shows  that  this  prediction  for  Z  is  quite  accurate  and 
thus,  one  would  expect  Eq.(25)  to  represent  closely  the  change  in  open 
area  resulting  from  linear  extension  of  the  threads.  Since  C  has  a 
value  on  the  order  of  10“^,  according  to  the  MIT  measurements, 

K2/K1  ~  10  which  means  that  AP  must  be  quite  large  (  ~  1000)  to 

produce  changes  in  area  of  10%.  For  most  practical  purposes,  then, 
one  may  neglect  this  effect. 

C.  CHANGE  IN  EFFECTIVE  OPEN  AREA  DUE  TO  ROTATION  OF 
THREADS  WITH  LOAD 

In  the  unstrained  condition  with  the  actual  pore  making  some 
angle  0  with  the  normal  to  the  cloth,  the  flow  must  turn  0°  on  the 
average  before  entering  the  pore.  Such  a  turn  will  cause  a  pressure 
drop  in  the  flow  so  that  the  force  driving  the  air  through  the  cloth  is 
not  Ap  but  rather  some  smaller  force.  Assuming  that  this  drop  in 
driving  force  does  not  depend  upon  Reynolds  number  and  that  the  approach 
'velocity  is  always  relatively  small,  it  appears  from  the  hydraulic  liter¬ 
ature  that  the  drop  in  effective  Pj  may  be  written 

APi  =  Ap  f(s)  sin  0  (27) 
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Inspection  of  data  in  Ref.  8  indicates  that  for  moderate  values 

of  Ap 

t- 

sin  0  ~  ^ (28) 

Apy 

where  y  is  0.2  to  0.5  or  so.  Thus,  on  the  average, 

APj^  «  fl(s)  AP*”^  (29) 

The  effective  change  in  open  areas  may  therefore  be  written 

1 

^  '  [  1  +  £i(s)AP-3]  (30) 

When  X  =  0.5,  Af  /A£  ~  AP®'^^  which  agrees  fairly  well  with  the 

trends  shown  by  Eq.  (18).  It  would  be  expected  that  the  losses  repre¬ 
sented  by  Eq.  (30)  will  be  the  only  additional  effect  of  pore  inclination 
since  it  appears  that  whatever  dynamic  head  exists  in  the  flow  issuing 
from  the  pores  is  lost  in  mixing  regardless  of  the  pore  inclination. 

The  analysis  given  above  is,  of  course,  of  a  very  qualitative 
nature  and  not  suitable  for  quantitative  determinations.  It  serves  to 
point  out,  however,  that  the  major  elastic  effect  in  parachute  cloth  so 
far  as  permeability  is  concerned  is  the  rotation  of  the  pores  under  load. 
From  measurements  available  of  this  rotation,  it  seems  reasonable  to 
conclude  that  the  elastic  contribution  to  x  at  low  AP  is  0.15  to  0.2  for 

^  4.  At  larger  AP  the 
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elastic  contribution  approaches  that  given  by  Eq.  (25),  A  more  complete 
analysis  of  the  effects  of  pore  rotation  was  not  possible  be  cans  e  of 
limitations  in  the  time  available  but  it  is  an  area  where  very  fruitful 
efforts  can  be  applied  in.  the  future.  Once  the  effects  of  pore  rotation 
are  established  quantitatively,  it  will  be  relatively  simple  matter  to 
determine  the  viscous  aspects  of  flow  through  parachute  cloth, 

D.  EFFECTS  OF  LOW  SOLIDITY  ON  FLOW  CHARACTERISTICS 

When  the  solidity  of  the  fabric  is  less  than  95%,  the  approach 
velocity  will  begin  to  be  significant.  This  has  the  effect  of  shifting  the 
critical  pressure  toward  1.0  as  indicated  in  Fig,  3  and  rh/A  will  be 
higher  for  the  same  AP.  The  reduction  in  solidity  also  permits  the 
character  of  the  approach  flow  to  change  from  incompressible  to 
compressible.  For  supersonic  approach  velocities,  the  effect  of 
adjacent  filaments  on  the  flow  over  any  given  filament  or  thread  can 
be  reduced.  This  happens  when  the  flow  between  the  threads  is  super¬ 
sonic.  The  maximum  solidity  for  which  this  is  possible  is  shown  in 
Fig.  4.  Under  these  conditions  the  drag  coefficient  of  the  screen  may 
be  considerably  less  than  with  subsonic  approach  velocities,  and  the 
permeabilities  considerably  higher  for  the  same  Because  the  shock 

detachment  distance  is  small,  however,  the  pressure  drop  across  the 
screen  measured  in  the  conventional  fashion  will  be  negative.  For  such 


IV.  EXPERIMENTAL  APPARATUS 


The  test  apparatus  is  shown  in  Figs.  5,  6,  7,  and  8.  The  figures 
are  for  the  most  part  self-explanatory.  Some  comments  on  the  design 
philosophy  and  experience  with  the  use  of  the  equipment  are  perhaps  in 
order  to  serve  as  food  for  thought  for  future  investigators. 

A.  MASS  FLOW  METERING 

Reference  to  the  results  will  show  that  for  these  tests,  mass 
flow  varies  ovor  a  range  of  10^.  The  tolerable  error  was  on 

the  order  of  2%  of  actual  indication  over  the  entire  range.  Originally, 
it  had  been  planned  to  use  two  venturi  meters  with  throat  diameters  of 
approximately  0.45”  and  0.9"  for  this  purpose.  ’  Venturi  meters  were 
favored  because 

1.  A  very  large  range  of  mass  flows  can  be  covered  accurately 
with  a  single  instrument,  provided  the  throat  velocity  remains  high,  i.e. 
the  pressure  drop  is  large  enough  to  be  read  accurately. 

2.  Very  little  pressure  drop  must  be  accepted  compared  to  an 
orifice  meter.  This  is  important  when  the  maximum  pressure  of  the 


system  is  limited. 


3.  The  reduction  of  data  is  simpler  than  with  an  orifice  meter 
because  the  flow  characteristics  are  better  understood. 

As  it  developed,  the  larger  venturi  meter  was  abandoned  because 
it  could  not  be  calibrated  easily  with  the  pumping  capacity  available  and 
because  the  smaller  meter  would  measure  the  maximum  flow  which  the 
nitrogen  supply  would  deliver  without  excessive  manometer  lengths. 
Considerable  time  and  effort  had  to  be  expended  in  shaping  the  remaining 
venturi  throat  to  produce  pressure  ratios  in  agreement  with  theory.  This 
was  necessary,  however,  to  take  advantage  of  the  simplicity  afforded  by 
the  theory.  To  avoid  disturbing  the  flow,  very  small  (.028")  pressure 
source  holes  were  needed.  To  keep  the  lag  within  reasonable  bounds, 
a  number  (4)  of  these  holes  were  drilled  around  the  periphery.  Because 
of  very  small  drifts  the  axial  locations  of  the  holes  during  drilling  the 
measured  pressure  at  each  hole  was  slightly  different  (+  2%).  By  using 
only  those  taps  which  averaged  nearest  0.528  Pj  for  choked  flow,  it  was 
possible  to  measure  the  choked  flow  value  within  the  reading  accuracy  of 
the  manometers. 

A  well-made  venturi  meter  should  require  calibration  only  at  low 
Reynolds  numbers.  This  simplifies  the  calibration  task  considerably 
since  the  capacity  of  the  calibrating  meter  or  device  need  not  be  large. 
This  assumption  seems  to  be  born  out  in  the  calibration  results  shown 

in  Fig.Al.  The  scatter  in  the  calibration  at  very  low  Reynolds  number 
is  undesirable,  but  since  the  device  was  not  used  at  these  conditions,  it 
is  not  important  to  the  data  reduction. 
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During  operation  with  fabrics  under  test,  it  was  found  that  the 
mass  flows  and  pressure  levels  were  such  that  the  venturi  rrseter  could 
not  be  used  accurately  for  mass  flows  less  than  0.5  gm/ sec,  except 
with  the  very  porous  stainless  steel  fabric.  Accordingly,  rotometers 
were  used  to  measure  m  down  to  0.01  gm/sec. 

The  rotometers  were  found  to  be  within  guaranteed  accuracy 
only  for  supply  pressures  above  atmospheric.  With  a  supply  pressure 
limit  of  30  psig,  each  rotometer  could  therefore  cover  about  a  decade 
in  mass  flow.  Thus,  two  rotometers  were  required  to  cover  the  range. 
The  first  covered  the  range  from  0.011  gm/sec  to  0.25  gm/sec.  The 
second  covered  the  range  from  0.1  gm/sec  to  2.4  gm/sec. 

The  rotometers  were  not  calibrated  further.  The  prediction 
charts  supplied  by  the  manufacturers  were  used  to  compute  curves  of 
indication  vs.  mass  rate  for  a  range  of  supply  pressures  and  temper¬ 
atures.  During  testing,  efforts  were  made  to  operate  on  these  curves. 
Individual  rotometers  are  guaranteed  by  the  manufacturer  to  be  within 
2%  of  these  curves.  This  was  considered  acceptable. 

Temperatures  used  in  computing  the  mass  flow  were  measured 
with  a  sensitive  ASTM  mercury  bulb  thermometer  inserted  into  the 


stream. 


B,  PRESSURE  MEASUREMENT 

Because  of  the  need  for  accuracy  over  a  wide  range,  at  least 
two  instruments  were  required  for  each  pressure  except  that  upstream 
of  the  rotometer.  This  was  measured  with  a  single  gauge  having  0.1  psi 
divisions.  The  venturi  meter  employed  a  mercury  U-tube  with  one  side 
open  to  atmosphere  (measured  by  a  mercury  barometer),  a  Wallace  and 
Tiernan  6"  dial  gauge  (0-100  mm,  l/2  mm  divisions)  accurate  to 
+  0.3  mm  and  a  Todd  Scientific  McLeod  gauge  (0-.5  mm,  0-5  mm, 

0-25  mm  scales)  to  measure  the  upstream  pressure.  The  differential 
pressure  was  measured  by  a  mercury  and  a  silicone  oil  (Dow  Corning 
DC-200,  sp.gr.  =  .97)  manometer.  The  silicone  oil  has  a  very  low 
vapor  pressure  but  is  hygroscopic;  frequent  cleaning  of  manometer  oil 
from  the  pressure  lines  was  necessary  because  of  violent  outgassing 
which  occurred  whenever  the  system  was  exposed  to  atmosphere  for 
any  length  of  time.  As  a  result,  the  system  was  later  kept  under  vacuum 
at  all  times  except  during  model  changes  and  repairs. 

The  test  sample  pressures  were  measured  in  a  similar  fashion. 
Provisions  were  made  on  both  panels  for  sealing  off  the  McLeod  gauges 
and  the  dial  gauges  for  tests  when  their  ranges  were  exceeded.  This 
proved  necessary  inspite  of  other  possible  arrangements  when  the  dial 
gauges  were  damaged  during  leak  checking,  and  mercury  from  the 
McLeod  gauge  was  blown  into  the  small  vacuum  pump.  A  trap  was  later 
provided  for  this  pump. 
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The  manometer  valving  arrangement  was  developed  when  other 
simpler  techniques  could  not  be  operated  with  s\ifficient  rapidity  or 
reliability  to  prevent  the  oil  manometers  from  blowing  over.  The 
valves  are  ordinary  bronze  steam  cocks.  They  proved  to  require  constant 
attention  to  keep  them  reasonably  leak  free,  however.  Lower  pressure 
tests  will  obviously  require  a  completely  different  type. 

Originally,  netjprene  rubber  tubing  was  used  for  the  pressure 
lines,  but  this  proved  too  susceptible  to  small  leaks  and  degassing  and 
it  was  replaced  by  copper  where  practical. 

Even  with  these  precautions,  it  developed  that  small  pressure 
differentials  (  <  2  mm),  particularly  at  low  ?£,  were  difficult  to  measure 
because  of  the  relatively  large  volume  of  the  measuring  system  and  the 
small  mass  flow  going  through  the  sample.  If  tests  are  ever  to  be  confined 
to  the  low  Ap,  low  P2  regime,  design  changes  are  indicated. 

It  should  be  noted  in  passing  that  in  designing  the  apparatus  and 
associated  instrumentation  for  this  experiment,  the  effort  was  made  to 
keep  it  as  simple  and,  hopefully  therefore,  as  reliable  and  inexpensive 
as  possible.  In  research  work,  some  changes  are  always  necessary  no 
matter  how  careful  the  initial  design.  It  was  felt,  therefore,  that  the 
procedure  adopted,  while  undoubtedly  requiring  more  changes  than  a 
more  elaborate  system,  was  in  the  long  run  just  as  accurate  and  reliable 


and  considerably  less  expensive. 


C.  TEST  SECTION 


As  elaborated  more  fully  in  Section  V,  the  test  section  was 
designed  to  test  as  large  a  sample  as  possible,  to  provide  a  more  uniform 
velocity  profile  than  a  straight  tube,  and  to  keep  the  flow  approximately 
normal  to  the  cloth,  i.e.,  prevent  spreading  downstream  of  the  sample. 
Originally,  a  2.5"  test  sample  was  used  but  because  the  mass  flow  did 
not  fall  off  as  rapidly  with  ?£  as  had  been  expected,  this  area  proved  to 
be  too  large  for  the  pumping  capacity  available  at  low  ^2'  -Accordingly, 
a  1.75"  diameter  test  section  was  constructed.  This  was  used  for 
most  of  the  tests.  For  the  very  porous  materials,  however,  even  this 
proved  too  large  to  get  large  pressure  drops  and  a  0.20"  test  section 
was  constructed.  To  avoid  excessive  length  in  the  recovery  section 
downstream,  a  sudden  expansion  to  1.75"  was  accepted.  This  entailed 
the  possibility  of  jet  spreading  and  results  different  from  those  with  the 
other  sections.  A  thin  rubber  sheet  was  used  to  confine  the  open  area  to 
0.20"  diameter  should  the  sample  be  blown  away  from  the  face  of  the 
upstream  section. 

These  test  sections  were  constructed  from  solid  aluminum.  For 
the  hot  tests,  the  test  section  was  constructed  of  copper  sheets  in  a 
manner  which  provided  space  for  circulation  of  a  coolant  (water).  Some 
cooling  was  necessary  to  protect  the  vacuum  seals  from  the  heat  and  to 


protect  the  operator  and  other  portions  of  the  equipment.  As  discussed 
later,  however,  the  heat  transfer  to  wall  was  too  large  and  a  significant 
temperature  gradient  existed  at  high  temperatures. 

For  tests  below  750°  F,  an  ASTM  mercury  bulb  thermometer 
placed  in  the  center  of  the  flow  upstream  of  the  contraction  was  used  for 
temperature  measurement.  For  tests  above  this  temperature,  a  copper 
constantan  thermocouple  was  used. 


D.  HEATER 


The  heater  consisted  of  48  coils  of  nichrome  wire  inside  quartz 
tubes.  These  tub  eSj  dimpled  at  the  upstream  end  to  prevent  motion 
downstream,  were  held  loosely  in  part  baffles.  This  arrangement 
required  almost  all  the  flow  to  pass  through  the  tubes.  The  walls  of 
the  flow  channel  were  made  of  thin  stainless  steel  insulated  from  the 
steel  case  by  MgO.  Changes  in  heat  output  were  effected  by  changing 
the  voltage  on  the  coils. 

The  heater  operation  was  generally  satisfactory  except  for  the 
long  lag  time  at  low  mass  flow.  At  mass  flows  less  than  1  gm/sec, 
cool  down  time  from  250°  F  was  on  the  order  of  2  hours.  Since  the 
heater  was  designed  to  meet  the  high  power  requirements  associated 

with  high  temperatures  and  large  mass  flows,  this  lag  is  unavoidable. 
Broken  tubes  and  wires  resulted  on  occasion  when  the  mass  flow  was 
changed  suddenly  during  operation  at  very  high  temperatures. 
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E.  COOLER 


The  cooler  was  a  conventional,  straight  gas-tube  type  and 
functioned  satisfactorily  under  all  heat  inputs  with  minimum  pressure 
drop. 


F.  PUMP 


The  laboratory's  130  CFM  vacuum  pump  was  used  for  most  of 
the  tests.  Near  the  end  of  the  program,  a  larger  350  CFM  pump  was 
placed  in  service.  The  apparatus  was  modified  to  make  use  of  this 


added  capacity  and  a  few  of  the  later  runs  utilized  both  pumps. 


V.  EFFECT  OF  APPARATUS  GEOMETRY  ON 
PERMEABILITY  MEASUREMENTS 


When  the  apparatus  for  the  present  investigation  was  designed  at 
the  beginning  of  the  program,  the  influence  of  the  apparatus  geometry  on 
the  test  results  was  not  fully  appreciated.  The  object  had  been  merely 
to  provide  a  reasonably  uniform  approach  velocity  without  excessive 
pressure  drops  and  to  maintain  the  flow  normal  to  the  cloth.  After  the 
calibration  of  the  venturi  meter,  however,  it  became  clear  that  the 
velocity  distribution  of  the  flow  as  it  approaches  the  test  sample  deter¬ 
mines  the  mass  flow  passed  by  the  cloth  for  a  given  AP  and  P2.  This 
can  be  seen  very  readily  from  the  following  considerations.  If  one 
assumes  that  the  test  channel  walls  are  parallel,  then  the  pressure 
normal  to  the  wall  is  constant  across  the  channel.  The  velocity  profile, 
however,  may  range  from  nearly  uniform  to  fully  parabolic.  When  a 
parabolic  profile  exists,  most  of  the  mass  flow  is  concentrated  near  the 
center.  Because  the  pressure  drop  across  the  cloth  increases  with 
increasing  mass  flow,  the  flow  near  the  center  will  have  a  higher  than 
average  pressure  drop.  Since  the  drop  must  be  the  same  across  the 
entire  cloth,  the  net  effect  is  that  for  a  parabolic  profile  the  pressure  drop 
will  be  higher  than  for  a  uniform  profile  with  the  same  mass  flow. 
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The  deficit  in  miass  flow  for  a  given  AP  and  P2  would  seem  to 
be  related  to  the  discharge  coefficient  of  the  test  channel.  Assuming  that 
the  test  channel  in  the  present  experiments  -  because  of  its  geometric 
similarity  to  the  venturi  meter  -  can  be^analyzed  by  the  method  described 
in  the  Appendix  and  that  B  is  approximately  the  same  as  for  the  venturi 
meter,  one  may  then  write  the  mass  flow  which  would  be  passed  for  a 
given  Ap  and  P2  if  the  profile  is  uniform  in  terms  of  the  actual  mass 
flow  as  follows: 

^uniform  ~  ^actual  3.55  Vpi  ^ ’^actual 
This  relation  is  shown  graphically  for  the  test  sections  used  in  these 
experiments  in  Pigs.  9.  Because  of  the  large  mass  flows  necessary  to 
obtain  reasonable  pressure  drops  along  the  channel,  it  was  not  possible 
to  verify  this  relation  directly. 

The  test  setup  used  by  the  MIT  group  would,  on  the  basis  of  its 
design,  also  be  expected  to  reduce  the  mass  flow  through  the  sample 
compared  with  that  which  would  be  passed  if  the  approach  velocity  were 
uniform.  Although  the  present  writer  has  less  theoretical  or  experi¬ 
mental  basis  for  an  accurate  estimate  of  the  characteristics  of  the  MIT 
apparatus  than  his  own,  he  is  inclined  to  believe  that  the  curves  in  Figs.  9 
provide  a  reasonable  first  approximation.  Many  other  permeability 
apparatus  probably  are  also  amenable  to  such  relatively  simple  analysis. 
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The  apparatus  used  in  the  tests  at  Georgia  Institute  of  Technology 
seems  to  be  in  a  different  category  with  respect  to  estimating  its  contri¬ 
bution  to  the  test  results.  In  this  apparatus,  the  metering  section  is 
followed  by  an  expanding  channel  which  is  terminated  by  the  test  sample. 
Although  the  test  sample  itself  is  relatively  large  (4"  diameter  in  one 
version  and  4"  square  in  another),  there  may  be  significant  nonuniformity 
in  the  flow  resulting  from  separations  from  the  wall  in  the  expanding 
section.  These  nonuniformities  are  virtually  impossible  to  estimate 
without  detailed  information  on  the  geometry,  turbulence  level,  and  mass 
flow.  Thus,  one  can  conclude  only  that  while  there  may  be  effects  due  to 
flow  nonuniformities  in  this  apparatus,  they  are  undetermined  at  present. 

The  discussion  thus  far  has  dealt  only  with  the  geometry  of  the 
channel  upstream  of  the  test  specimen.  It  can  be  shown,  however,  that 
the  geometry  downstream  of  the  test  specimen  and  the  location  or  char¬ 
acteristics  of  the  gas  pump  also  contribute  to  the  results.  If,  for  example, 
the  channel  immediately  downstream  of  the  sample  enlarges  discon- 
tinuously,  or  even  very  rapidly,  the  downstream  pressure  may  be 
considerably  lower  for  the  same  mass  flow  than  if  the  diameter  remains 
constant.  Whether  this  happens  will  depend  on  the  dispersive  power  of 
the  sample  and  the  admittance  of  the  downstream  mass  sink  (atmosphere 
or  pump).  If  the  sample  produces  little  dispersion  in  the  stream,  or  if 
the  downstream  pressure  is  held  constant,  no  matter  what  the  mass  flow 
through  the  sample,  _toen  the  m  ~  AP  relation  will  be  independent  of 
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downstream  geometry.  On  the  other  hand,  if  the  sample  produces  a 
large  dispersion,  the  downstream  pressure  will  be  lower  for  the  same 
mass  flow  in  order  to  satisfy  mass  conservation.  This  drop  may  not 
materialize,  however,  if  it  will  result  in  changes  in  the  inlet  conditions 
of  a  constant  volume  pump.  Rather,  the  upstream  pressure  will  tend  to 
rise. 

It  must  now  be  evident  to  the  reader  that  permeabilities  of  given 
cloth  samples  in  the  present  apparatus  will  be  quite  different  from  those 
measured  in  the  MIT  apparatus.  Generally  because  the  approach  velocity 
is  probably  more  nonuniform  and  the  channel  enlarges  discontinuously 
downstream  of  the  test  sample,  the  MIT  data  would  be  expected  to  show 
a  larger  pressure  drop  for  the  same  mass  flow.  The  MIT  data  can  also 
be  expected  to  show  a  larger  percentage  drop  in  mass  flow  with  decrease 
in  P2  than  the  present  data. 

While  one  may  argue  that  flow  through  a  parachute  canopy  i  s 
certainly  dispersed,  it  woxild  -seem  that  cloth  permeability  measurements 
per  se  should  be  carried  out  in  a  constant  area  stream  tube  w’ith  a  uniform 
approach  velocity  which  is  normcd  to  the  sample  surface.  Study  of 
pressure  distributions  on  complete  parachute  canopies  reported  in  Ref.  11 
suggest  that  this  is  the  most  reasonable  basis  for  generalizing  test  data 
on  small  cloth  samples  to  full  scale  canopies. 
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VI.  EXPERIMENTAL  PROCEDURE 
AND  HANDLING  OF  DATA 


A.  COLD  FLOW  TEST  PROCEDURE 

Runs  were  usually  begun  at  a  downstream  pressure  of  approxi¬ 
mately  760  mm  Hg.  The  downstream  fine  metering  valve  was  opened 
slightly  and  dry  nitrogen  was  admitted  through  one  of  the  rotometers. 
The  pressure  upstream  of  the  rotometer  was  adjusted  commensurate 
with  the  upstream  temperature  so  that  the  density  through  the  rotometer 
was  1,  2,  or  3  atmospheres.  The  pressure  drop  across  the  cloth  was 
recorded.  Oil  manometers  were  used  for  Ap  <  60  mm  Hg.  By 

opening  the  downstream  metering  valve  and  increasing  the  gas  in-flow 
rate,  the  downstream  pressure  was  held  at  760  mm  and  the  AP 
increased.  At  mass  flows  higher  than  1  to  2.5  gm/sec,  the  flow  was 
diverted  from  the  rotometers  and  passed  directly  through  the  venturi. 
Upstream  venturi  pressure  was  read  with  a  mercury  manometer  for 
P^  >  100  mm  Hga  and  with  a  Wallace  and  Tiernan  dial  gauge  for 

Pj  <  100  mm  Hga.  AP  was  read  with  the  silicone  oil  manometer 

for  Ap  <  60  mm  and  with  the  mercury  manometer  for  AP  >  60  mm. 

Points  were  taken  more  or  less  at  random,  the  effort  being  to 
take  sufficient  points  to  permit  fairing  through  the  usual  random  vari¬ 
ations  and  to  keep  the  log  AP  interval  fairly  constant.  The  upper  limit 
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0 


■was  determined  by  which  came  first:  (a)  gas  supply  limitations  of 

65  gm/lec  max,  (b)  pump  capacity  limitation  of  approximately 

\ 

0.1  P2  mm  Hg  gm/sec,  (c)  manometer  deflection  limitations  of 
990  mm  Hg,  or  (d)  fabric  failure. 

Following  the  completion  of  a  run  at  =  760  mm,  the  down¬ 
stream  pressure  was  decreased  and  the  procedure  repeated,  usually 
for  descending  AP.  The  third  P2  run  utilized  an  ascending  AP,  etc. 
The  selection  of  P2  was  made  to  facilitate  definition  of  the  trends  with 
P2.  At  the  low  P2's,  the  McLeod  gauge  was  used  to  check  the  pressures. 

Usually  a  second  or  a  third  sample  was  spot  checked  to  determine 
the  degree  of  scatter  from  the  results  of  the  first  sample.  Four  different 
test  apparatus  were  available;  a  particular  material,  however,  was  tested 
in  only  three.  For  the  most  part,  different  samples  were  used  in  the  tests 
in  the  different  test  sections. 

Perhaps  it  might  be  well  to  comment  at  this  point  on  the  use  of 
dry  nitrogen  as  the  test  gas  in  place  of  air.  First  of  all,  nitrogen  was 
available  in  the  laboratory  and  a  dry  air  supply  suitable  for  mass  flows 
less  than  30  gm/sec  but  more  than  1  gm/sec  was  not.  Secondly,  nitrogen 
is  very  similar  to  air  in  density  and  viscosity  so  that  the  results  should 
be  comparable.  Thirdly,  the  life  of  open  wound  and  other  simple  heating 
elements  is  much  longer  in  nitrogen  than  in  air  which  made  the  design 
and  construction  of  the  heater  easier  and  less  costly.  The  only  area 
where  the  use  of  nitrogen  may  make  a  difference  is  in  the  high 
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temperature  degradation  of  organic  fibers.  There  is  evidence  in  the 
literature  to  indicate  that  oxygen  causes  more  rapid  degradation  than 
nitrogen.  This  point  was  left  for  future  investigation,  however. 


B.  HOT  FLOW  PROCEDURE 

Because  of  the  difficulty  in  regulating  the  output  of  the  heater, 
the  data  at  a  given  temperature  level  were  taken  at  constant  mass  flow. 

The  variation  in  pressure  drop  with  changes  in  P2  was  recorded.  The 
temperature  was  then  increased  and  when  a  stable  value  was  reached, 
the  procedure  repeated.  The  temperature  was  increased  in  increments 
until  failure  was  observed  or  to  a  planned  1500°  F,  whichever  occurred 
first.  The  temperature  increments  were  chosen  in  an  effort  to  define 
the  effects  of  temperature  with  the  minimum  number  of  points. 

C.  DATA  REDUCTION  PROCEDURE 

Except  for  the  viscous  corrections  applied  to  the  measured  mass 
flow  and  the  correction  for  approach  velocity  nonuniformity,  the  basic  data 
reduction  procedure  is  self-indicating.  The  viscous  corrections  applied 
to  the  venturi  mass  flow  were  obtained  from  the  venturi  calibration  (see 
Appendix  and  Fig.  Al).  The  rotometer  data  was  used  directly.  All  the 
results  so  obtained  are  shown  in  Figs.  10  through  ZO.  Later,  the 
estimated  effects  of  approach  velocity  nonuniformity  were  applied  to  the 
faired  curves  and  these  were  then  redrawn  as  Figs.  21  through  27. 


VII.  DISCUSSION  OF  RESULTS 


A.  COLD  FLOW  TESTS 

s 

1 .  Variations  from  Sample  to  Sample 

Generally,  the  thinner  and  more  rigid  the  sample,  the  less  the 
variation  in  permeability  from  sample  to  sample.  The  periiieability  of 
the  stainless  steel  fabric  was  virtually  the  same  for  all  samples.  Devi¬ 
ations  were  usually  confined  to;  the  scatter  band  for  any  particular  run. 
The  Mil-C-7350B  Type  I  fabric  also  gave  reasonably  good  agreement 
from  sample  to  sample.  The  MU-CtSOZIB  Type  II  fabric  showed  vari¬ 
ations  on  the  order  of  20%  while  the  Mil-F-9084  Type  VIII  fabric  showed 
variations  of  around  50%  from  sample  to  sample.  The  HT-1  fabric 
showed  the  same  variations  of  the  Mil-C-7350B  Type  I  fabric. 

It  is  important  to  note,  however,  that  the  ratio  of  permeability 
at  one  downstream  pressure  to  that  of  another  for  the  same  AP  was 
nearly  the  same  for  all  samples  of  the  same  material,  indicating  that 
while  there  may  be  minor  differences  in  construction  of  the  various 
samples,  these  do  not  affect  the  qualitative  air  flow  characteristics. 

The  resiilts  of  other  tests  on  the  same  fabric  can  thus  be  converted  to 
higher  altitude  results  by  multiplying  the  present  data  by  the  ratio 

^760  test 
^760  present  data 
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2.  Effect  of  Sample  Area 

The  two  heavier  materials  were  tested  in  both  the  2.5"  and 
1.75  "test  sections.  After  correcting  the  results  for  flow  nonuniformity, 
it  was  found  that  the  mass  flow  per  unit  area  was  the  same  within  20% 
for  the  8021B  fabric  and  50%  for  the  9084  fabric.  This  would  indicate 
that  (a)  the  smaller  sample  size  is  sufficient  for  good  average  results, 
and  (b)  that  the  corrections  for  flow  nonuniformity  were  proper. 

3.  Comparison  with  Results  of  Other  Investigations 

Figs.  28,  29,  30,  and  31  give  comparisons  between  the  uncorrected 
(i.  e.  ,  for  approach  velocity  nonuniformity  and  downstream  spreading) 
results  of  several  investigations  of  the  same  material.  At  sea  level,  the 
use  and  MIT  tests  of  the  7350B  Type  I  fabric  show  approximately  the 
same  trend,  but  a  different  "gain",  the  difference  being  somewhat  more 
than  that  found  between  samples  of  7350B  Type  I.  It  seems  reasonable 
to  expect  that  this  is  a  result  of  downstream  spreading  in  the  MIT  appa¬ 
ratus.  It  will  be  noted  that  the  Georgia  Institute  of  Technology  apparatus 
gave  results  which  agreed  more  closely  with  the  USC  results  in  magnitude, 
but  with  a  slope  indicative  of  either  different  elastic  properties  or  a  change 
in  the  velocity  distribution  approaching  the  cloth. 

At  the  higher  altitudes  the  results  obtained  with  the  MIT  apparatus 
show  a  much  greater  drop  in  G  and  a  higher  value  of  x  than  with  the 
USC  apparatus.  This  is  to  be  expected  since  the  approach  velocity  may 


be  somewhat  more  nonutiiform  in  the  MIT  apparatus  and  the  downstream 
spreading  can  be  expected  to  increase  with  decreasing  Reynolds 
number. 

The  sea  level  results  obtained  with  the  8021B  fabric  are  somewhat 
anomalous.  The  gain  of  USC  results  is  only  about  70%  of  the  MIT  results, 
although  the  slopes  are  the  same.  From  tests  of  other  samples  of  the 
8021B  fabric  it  would  be  expected  that  the  curve  marked  USC  represents 
the  lower  limit.  Even  if  one  adds  20%  or  30%  the  results  will  still  fall 
25%  below  the  MIT  data.  The  material,  being  quite  dense,  would  tend  to 
promote  downstream  spreading  and  thus,  the  results  are  even  less  under¬ 
standable.  The  GIT  results  because  of  their  peculiar  slope  must  be 
regarded  as  affected  by  increasing  separation  from  the  walls  of  the 
stream  approaching  the  test  sample. 

At  the  higher  altitudes,  the  effect  of  apparatus  geometry  on  the 
MIT  results  as  opposed  to  the  USC  results  is  evident. 

4.  Predictability  of  Results 

Fig.  3  shows  the  predicted  permeability  of  the  stainless  steel 
fabric  based  on  the  two  cylinder  flow  analogy.  Comparison  with  the 
measured  results  (Figs.  11  and  21)  yields  excellent  agreement  for 
P2  >  46  mm  (about  63,000).  For  P^  —  46  the  measured  data  is 

increasingly  higher  than  the  predicted,  indicating  that  the  predicted 
viscous  effects  did  not  materialize.  At  P2  =  46mm,  the  Knuds en 


number  of  the  fabric  is  about  0.03  and  thus,  one  would  expect  that  the 
characteristics  will  begin  to  deviate  from  those  of  continuum  theory. 
However,  measurements  of  the  characteristics  of  cylinders  in  the 
transition  region,  between  continuum  flow  and  free  molecule  flow,  are 
not  yet  well  enough  established  to  permit  inclusion  of  these  effects  in  the 
prediction.  The  present  predictability  limit  is  thus  about  50,000'.  For 
altitudes  below  this,  it  has  been  demonstrated  that  the  two  cylinder  model 
gives  excellent  results  for  solidities,  at  least  as  high  as  .52.  Morgan's 
results  (Ref.  20)  show  the  same  feature  at  higher  Reynolds  numbers. 

Fig.  3  was  also  compared  with  the  measured  permeabilities  of  the 
other  materials.  The  qualitative  agreement  with  the  7350B  Type  I  fabric 
was  good  for  P2  ^  46  mm,  but  some  elastic  effects  or  altered  viscous 

effects  did  produce  some  differences.  It  is  difficult,  as  noted  earlier, 
to  assign  the  cause  to  one  or  the  other.  The  trends,  however,  seem  in 
closer  agreement  with  the  two-cylinder  model  than  with  either  the  smooth 
nozzle  or  the  orifice. 

The  permeabilities  of  the  two  heavier  fabrics  showed  somewhat 
poorer  agreement  with  the  qualitative  predictions  for  the  screen  at 
P2  >  46  mm,  indicating  that  the  viscous  effects  may  be  somewhat 
different  than  assumed. 

For  P2  <  46  mm,  the  permeabilities  of  all  fabrics  evidenced 

rarefaction  effects.  It  appears,  therefore,  that  satisfactory  prediction 
of  the  permeabilities  of  fabrics  at  high  altitudes  must  await  a  satisfactory 
description  of  their  transition  flow  characteristics. 


45 


5.  Empirical  Description  of  Change  in  Permeability  with  Altitude 

An  attempt  was  made  to  obtain  from  the  experimental  results 
the  functional  dependence  of  G  and  x  on  AP,  and  m,  but  no 

satisfactory  correlation  has  been  obtained  so  far.  For  the  present,  one 
must  accept  the  corrected  experimental  results  in  lieu  of  a  valid  scaling 
procedure  or  a  complete  analytical  description  of  the  permeability. 


B.  HOT  FLOW  TESTS 


At  temperatures  below  the  melting  point  of  the  material,  one 
would  expect  that  either  the  material  will  stretch  slightly  under  load 
or  that  it  will  begin  to  yield,  in  which  case  it  is  impossible  to  obtain 
steady  state  permeability  indications.  Since  it  has  been  shown  that 
simple  extension  of  the  fibers  results  in  an  almost  negligible  increase 
in  permeability,  one  would  expect  the  mass  flow  at  a  given  AP  and 
P2  to  depend  almost  entirely  on  density.  For  experimental  convenience 
however,  the  tests  were  designed  to  measure  the  variation  in  AP  with 
T  at  constant  m  and  P^.  For  this  case  one  would  expect  that  since 

m  _  G' 

«  =  AP^  at  constant  P7, 


one  would  find  that 
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The  results  showed  the  exponent  of  Tg  to  be  somewhat  more  than 
predicted  in  every  case.  The  reason  for  this  is  that  the  heat  transfer 
to  the  test  section  wall  tends  to  make  the  mass  flow  distribution  over 


the  face  of  the  sample  more  uniform,  thus  making  the  value  of  x  slightly 
smaller.  That  this  is  the  case  seems  evident  from  an  examination  of  a 
typical  sample  which  has  been  exposed  to  a  temperature  much  higher 
than  its  melting  point  (Fig.  32).  It  will  be  seen  that  only  the  center  1  ” 
is  melted,  the  outer  3/8  of  an  inch  on  all  sides  is  unmelted. 

While  such  high  heat  transfer  is  undesirable  from  the  point  of 
view  of  data  interpretation,  the  analysis  above  and  the  experimental 
results  seem  to  demonstrate  conclusively  that  the  long  term  effects  on 
permeability  of  changes  in  stagnation  temperature  can  be  accounted  for 
entirely  by  multiplying  by  ^/Tg^/Tg^. 

\  The  8021B  fabric  was  tested  to  a  temperature  of  400°  F  without 
failure\  Tests  at  intermediate  values  of  275°  F  and  325°F  were  also 
made.  Temperature  overshoot  of  the  heater,  however,  made  it  impos¬ 
sible  to  approach  the  melting  temperature  more  closely  without  sample 
failure. 

The  new  organic  fiber  HT-1  was  tested  to  650°F  without  failure. 
Tests  at  intermediate  temperatures  of  400°  F,  500°  F  and  600° F  were 
also  run.  The  sample  failed  when  an  attempt  was  made  to  test  it  at 
700°  F.  The  sample  is  shown  in  Fig.  33. 


The  glass  fabric  was  tested  at  275°  F,  500°  F,  and  750°  F 
without  evidence  of  any  damage.  It  had  been  hoped  to  test  it  at  higher 
temperatures  but  failure  of  the  heater  in  the  last  days  of  the  testing 
program  prevented  this. 

The  stainless  steel  fabric  was  tested  at  750°  F  and  an  indicated 
1200°F.  There  is  reason  to  believe,  however,  that  the  thermocouple 
indication  may  have  been  too  low  (due  to  heat  losses  along  the  wires  and 
reradiation  to  the  wall).  This  temperature  was  the  maximum  on  the 
galvanometer  scale.  It  had  been  planned  to  use  another  indicator  for 
higher  temperatures,  but  failure  of  the  heater  at  this  point  forced 
abandonment  of  this  plan  since  repair  of  the  heater  would  have  required 
about  a  week  and  the  time  was  not  available. 

The  results  of  the  various  hot  tests  are  shown  in  Figs.  34  through 
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VIII.  CONCLUSIONS 


The  present  program  has  shown  for  the  first  time  that: 

1.  The  internal  geometry  of  the  permeometer  can  have  an 
important  effect  on  the  measured  permeability.  Results  most 
similar  to  those  for  a  portion  of  a  parachute  canopy  should  be 
obtained  with  a  permeometer  employing  a  bell-mouth  inlet,  with 
the  sample  placed  across  the  end  of  the  inlet  followed  by  a  tube 
having  the  same  internal  diameter  as  the  end  of  the  inlet. 

2.  Viscous  effects  are  present  under  all,  except  the  very 
highest,  mass  flows  for  all  materials  tested. 

3.  Rarefaction  effects  appear  to  be  present  at  altitudes 
above  60, 000' . 

4.  The  general  equation  for  mass  flow  through  the  cloth  is 


The  expression  f(P2/Plj  AP,  Re)  is  in  general  different 
for  each  type  of  material  and  cannot  as  yet  be  determined 
analytically.  For  constructions  where  the  dependence  on 
P2/P1  and  Ap  is  small,  satisfactory  predictions  can  be 
obtained  for  f(s,  Rg)  if  s  <  0.5  or  so  and  the  Knudsen 
number  is  less  than  .01. 
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5.  The  principal  change  in  permeability  resiilt^ng  from 
material  elasticity  is  connected  with  the  rotation  of  the  fabric 


pQ3-0g  \xnder  load*  Changes  associated  with  linear  extension 
of  the  yarn  arises  only  at  very  high  AP's. 

For  the  materials  tested,  the  effects  of  elevated  temperatures 
on  permeability  are  derived  entirely  from  changes  in  air  density. 


A.  RECOMMENDATIONS 


» 

1.  From  the  point  of  view  of  basic  research,  the  flow  charac¬ 
teristics  through  typical  cloth  samples  should  be  determined  under 
very  low  density  conditions  to  delineate  the  rarefaction  effects.  Tests 
should  be  conducted  with  P2  as  low  as  30^i  Hg  and  with  AP/P2  up  to 
about  4. 

2.  This  will  require  that  the  test  sample  be  at  least  4"  in 
diameter  to  avoid  excessive  approach  velocity  nonuniformity. 

3.  A  qualitative  analysis  -  both  theoretical  and  experimental  - 
of  the  effect  of  pore  inclination  on  its  mass  transfer  characteristics 
should  be  initiated.  This  should  be  coupled  with  an  investigation  to 
determine  the  discharge  characteristics  of  smooth  nozzles,  orifices 
and  two-cylinder  nozzles  at  low  densities. 

4.  Analysis  of  the  data  already  on  hand  should  be  continued  in 
an  effort  to  define  the  viscous  and  rarefaction  effects  more  satis¬ 
factorily. 

5.  Except  when  elastic  effects  are  directly  and  reversibly 
related  to  gas  temperature,  tests  at  elevated  temperatures  need  not 
be  conducted. 


51 


IX.  REFERENCES 


1.  Cornell,  W.  G.  ,  "Losses  in  Flow  Normal  to  Plane  Screens", 
Transactions  of  the  ASME  80,  p.  791,  1958. 

2.  Morgan,  P.  G.  ,  '.'High  Speed  Flow  Through  Perforated  Plates", 

Journal  of  the  Royal  Aeronautical  Society  64,  Feb.  I960. 

3.  Schlichting,  H.  ,  Boundary  Layer  Theory,  4th  ed.  ,  McGraw-Hill 
Book  Co.  ,  I960. 

4.  Lavier,  H.  W.  S. ,  "Air  Permeability  of  Parachute  Cloth",  WADC 
TR  52-283,  part  II,  August  1953. 

5.  Scheidegger,  Adrian  E.,  The  Physics  of  Flow  Through  Porous  Media. 
University  of  Toronto  Press,  1957. 

6.  Green,  Leon  andDuwez,  Pol,  "Fluid  Flow  Through  Porous  Media", 
J.  of  Appl.  Mech.  18,  No.  1,  pp.  39-45,  March  1951. 

7.  Lavier,  H.  W.  S.  and  Boteler,  W.  C.  ,  "Air  Permeability  of  Parachute 
Cloth",  WADC  TR  52-283,  part  V,  April  1955. 

8.  Seshadri,  Brown,  Backer,  Krizik,  and  Mellen,  "Air  Flow  Charac¬ 
teristics  of  Parachute  Fabrics  at  Simulated  High  Altitudes",  WADC 
TR  59-374,  August  1959. 

9.  Klein,  Lermond  and  Platt,  "Research  Program  for  the  Development 
of  a  Design  Procedure  to  Engineer  Parachute  Fabrics",  WADC 

TR  58-65,  May  1958. 

0.  Bickford,  Kuehl,  and  Rusk,  "Study  of  the  Control  of  Permeability  of 
Nylon  Parachute  Cloth  at  High  and  Low  Differential  Pressures", 

WADC  TR  54-468,  March  1955. 

1.  Heinrich,  Ballinger  and  Ryan,  "Pressure  Distribution  in  Transonic 
Flow  of  Ribbon  and  Guide  Surface  Parachute  Models",  WADC  TN 
59-32,  February  1959. 


12. 


Crocco,  Luigi,  "One-Dimensional  Treatment  of  Steady  Gas  Dynamics", 
Fundamentals  of  Gas  Dynamics,  H.  W.  Emmons,  editor,  Princeton 
University  Press,  1958. 

13.  Pai,  S.  I.,  Viscous  Flow  Theory,  Vol.  II,  D.  Van  Nostrand,  1957. 

14.  A.  S.  M.  E.  ,  Fluid  Meters. 

15.  Liepmann,  H.  W.,  "A  Study  of  Effusive  Flow",  Aeronautics  and 
Astronautics.  Pergarnon  Press,  I960. 

16.  Stearns,  Jackson,  Johnson,  and  Larson,  Flow  Measurement  with 
Orifice  Meters.  Van  Nostrand,  1951. 

17.  Liepmann,  Hans  W.,  "Gaskinetics  and  Gasdynamics  of  Orifice  Flow", 
Journal  of  Fluid  Mechanics,  Vol.  10,  part  1,  1961. 

18.  Perry,  J.  A.  ,  Jr.  ,  "Critical  Flow  through  Sharp-Edged  Orifices", 
Transactions  of  the  ASME,  71,  1949. 

19.  Marks  Mechanical  Engineers  Handbook,  McGraw-Hill  Book  Co.  , 

New  York. 

20.  Morgan,  P.  G.  ,  "High  Speed  Flow  through  Wire  Gauges",  Journal  of 
the  Royal  Aeronautical  Society,  Vol.  63,  August  1959. 

21.  Rouse,  Hunter,  Fluid  Mechanics  for  Hydravilics  Engineers,  Dover, 

1961. 

22.  Rouse,  Hunter,  Elementary  Mechanics  of  Fluids,  John  Wiley,  1947. 

23.  Gowen,  Forrest  E.  and  Perkins,  Edward  W.  ,  "Drag  of  Circular 
Cylinders  for  a  Wide  Range  of  Reynolds  Numbers  and  Mach  Numbers", 
NACA  TN  2960. 

24.  Eshbach,  editor,  Handbook  of  Engineering  Fundamentals.  John 
Wiley,  1952. 

25.  Kennard,  E.  H.  ,  Kinetic  Theory  of  Gases,  McGraw-Hill  Book  Co.  , 
1938. 


APPENDIX  I 


AN  EXTENDED  THEORY  OF  GAS  FLOW  RATE 
THROUGH  A  NOZZLE 

A.  GENERAL 

It  is  usually  assumed  that  the  flow  through  a  convergent-divergent 
nozzle  can  be  represented  quite  adequately  by  treating  the  fluid  as  if  it 
were  adiabatic  and  inviscid  with  uniform  properties  at  each  axial  location. 
It  has  been  shown  —  by  experiment  and  also  by  more  inclusive  analysis  — 
that  these  assumptions  are  well-substantiated,  provided  the  throat 
entrance  and  exit  geometries  are  not  too  abrupt,  the  Reynolds  number 
based  on  throat  diameter  is  not  too  low,  and  the  throat  Knuds en  number 
is  not  too  high.  The  equations  derived  under  these  assumptions  are 
well-known  and  therefore  will  not  be  repeated  here.  For  the  purposes 
of  the  present  investigation,  only  the  equation  relating  the  mass  flow 
and  the  pressures  at  the  throat  and  some  upstream  location  is  of  interest. 
It  is  the  purpose  of  the  following  analysis  to  determine,  if  possible,  the 
manner  in  which  this  relation  must  be  modified  to  treat  cases  for  which 
some  of  the  foregoing  assumptions  are  violated. 
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B.  MODIFICATION  FOR  REDUCED  REYNOLDS  NUMBER 


Since  gases  are  in  fact  slightly  viscous,  reducing  the  Reynolds 
number  can  be  expected  to  produce  increasing  deviations  from  the  flo-w 
rate  predicted  by  the  ideal  theory.  In  contrast  to  the  widespread  know¬ 
ledge  of  the  ideal  theory  result,  there  appears  to  be  no  exact  solution 
for  Re^j  <<00  .  Crocco  (Ref.  12)  examined  the  situation  in  some  detail 

but  in  a  rather  approximate  fashion  which,  even  so,  does  not  permit 
simple  calculation  of  the  ratio  of  actual  mass  flow  rate  to  the  ideal  value, 
i.e.,  the  discharge  coefficient.  His  analysis  shows,  however,  that  the 
value  of  C(j  is  determined  almost  entirely  by  events  occurring  upstream 
of  or  at  the  throat.  The  exception  occurs  when  the  pressure  gradient 
downstream  of  the  throat  is  unfavorable  enough  to  produce  separation 
upstream  of  the  throat.  This  situation  is  discussed  in  detail  later. 
Schlichting  (Ref.  3)  and  Pai  (Ref.  13)  discuss  theoretical  and  experimental 
results  for  incompressible  flow  in  a  straight  pipe  under  both  laminar 
and  turbulent  flow  conditions,  and  both  also  give  the  exact  solution  for 
incompressible  viscous  flow  in  a  converging  two-dimensional  channel. 
Unfortunately,  while  such  cases  may  be  similar  in  many  ways  to  the 
present  one,  they  are  not  sufficiently  close  to  warrant  inferring  from 
them  the  quantitative  characteristics  of  compressible  nozzle  flow. 
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« 

In  the  absence  of  an  adequate  theoretical  treatment  one  can  consult 
the  hydraulics  literature  for  measurements  of  the  discharge  character¬ 
istics  of  venturi  meters,  nozzles  and  orifices.  These  are  usually  limited 
to  incompressible  flow,  however.  One  does  find  (Ref.  14)  intimations  that 
the  discharge  characteristics  are  the  same  for  a  compressible  fluid, 
provided  the  normalization  is  carried  out  using  the  inviscid  compressible 
flow  theory.  However,  no  data  was  presented  to  support  this  conjecture. 
Also,  recent  editions  of  the  standard  engineering  handbooks  have  tended 
to  regard  the  earlier  measurements  as  questionable,  but  because  of 
apparent  lack  of  interest  in  the  nozzle  as  a  standard  flow  rate  measuring 
device,  particularly  at  low  Reynolds  numbers,  have  included  little  new 
data.  They  suggest  individual  calibration  of  the  device  used  as  the  most 
satisfactory  approach. 

In  view  of  this  situation,  it  was  felt  that  a  detailed  calibration  of  the 
venturi  meter  used  in  this  apparatus  for  metering  the  higher  mass  flow 
would  provide  some  instructive  information  on  the  behavior  of  compres¬ 
sible  flow  in  nozzles  at  low  Reynolds  numbers  and  serve  as  a  standard 
against  which  to  compare  the  results  of  any  theory  which  might  be  proposed 
to  explain  it.  The  results  are  shown  in  Fig.Al.  The  fact  that  no  measur¬ 
able  change  in  vs.  Re^J  was  observed  when  Pg/Pi  >  0.528  seems  to 
confirm  earlier  beliefs  in  this  regard.  This  might  be  expected  since  the 
maximum  Mach  number  in  the  converging  section  is  1.0  and  the  displace¬ 
ment  thickness  on  a  flat  plate  at  least  varies  but  little  between  M  =  0 
and  M  =  1.0. 
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After  some  study,  it  was  found  that  the  experimental  results 
could  be  explained  satisfactorily  by  the  following  semi -empirical 
analysis: 

Consider  the  convergent  portion  of  an  axisymmetric  nozzle  as 
shown  in  the  sketch 


When  10^  <  Re^j  <  10^>  one  is  probably  justified  in  assuming  that  the 
flow  at  the  throat  consists  of  an  isentropic  core  and  a  laminar  boundary 
layer  with  displacement  thickness  S  The  discharge  coefficient, 
can  be  represented,  therefore,  by 

M  -  \  ^  48*  4  8*^ 

Cd  d  /  "  ^  ■  d 


When  2  8*/d  <<  1, 


1  + 


48^ 


(A2) 


If  one  assumes  that  8*  =  ^1  ^  ,  Eq.  (A2)  becomes 

Rex 


Cd 


i  +  15iiL 

d^Rex 


(A3) 
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seems  reasonable  since  one  finds 
this  Reynolds  number  dependence  for  both  a  flat  plate  and  a  circular 
cylinder  and  since  it  has  been  shown  that  two-dimensional  boundary 


layer  theory  gives  good  results  for  flow  inside  ducts  if 


<  <  1. 


For  a  given  nozzle,  x/d  =  constant,  K2.  |^It  should  be  noted, 
however,  that  if  the  contraction  ratio  is  >5  (in  terms  of  diameters), 

X  should  probably  be  measured  from  the  position  at  which  the  velocity 
first  becomes  significant,  i.  e.  ,  when  D/d  is  first  less  than  5.J  Using 
the  fact  that  K2  gives  the  relation  between  s  and  d,  one  may  write 


_s: 

d 


Kly^ 

Hence,  Eq.  (A3)  can  be  written 


(A4) 
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1  + 


y^ed 


or 


C.  = 


1  + 


B 


(A5) 


y^ed 


The  value  of  Kj,  and  hence  of  B,  depends  upon  dP/dx.  One  may  estimate 
the  value  of  from  the  exact  analysis  of  flow  over  a  cylinder.  In  that 
analysis 


8*  = 


0.8R 
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(A6) 
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at  the  90°  point.  For  R/d  =  1.66  »  K21  B  =  4.12.  Tests  of  a  nozzle 
■with  this  entrance  geometry  gave  B  =  4.0.  The  calibration  data  are 
shown  in  Fig.  (Al)'.  Dry  nitrogen  was  the  test  fluid. 

Eq.  (A5)  can  be  expected  to  yield  reasonable  agreement  with  exper¬ 
iment  down  to  Reynolds  number  for  which  the  boundary  layer  completely 
fills  the  throat.  Continuing  the  use  of  the  analysis  of  flow  over  a  cylinder, 
one  has  S*/8  =  0.3  at  the  throat.  This  result  assumes  a  two-dimensional 
flow,  however.  In  the  three-dimensional  case  of  a  nozzle,  the  pressure 
gradient  is  somewhat  more  favorable  and,  according  to  the  Karman- 
Pohlhausen  analysis,  S*/8  is  therefore  less  than  0.3.  If  one  takes 
SVS  =  0.27  then  one  has  =  0.65  which  agrees  with  experiment. 

As  the  Reynolds  number  is  decreased  below  this  point,  the  flow 
velocity  at  the  throat  decreases.  This  comes  about  because  the  pressure 
at  the  throat  does  not  decrease  -  as  near  as  can  be  determined  experi¬ 
mentally  -  compared  with  the  upstream  pressure  while  the  viscous 
dissipation  continues  to  increase.  Because  the  velocity  decreases  and 
the  dissipation  results  in  increased  temperature,  the  entire  flow  soon 
becomes  nearly  isothermal.  For  such  flows,  the  discharge  is  the  same 
as  for  incompressible,  viscous  flow  for  which  ih  ~  AP  in  a  straight 
pipe.  Isentropic  nozzle  theory,  however,  gives  m  ~ v^AP  for  constant 
P2  so  that  in  a  slowly  contracting  channel  na/iiiinviscid  ~  R-e'  Thus, 

Cd  = 


C 


(A7) 


The  value  of  C  depends,  of  course  on  dP/dx.  If  one  assumes  that  for 


a  unit  length  of  straight  pipe  AP  =  .472Pj^,  then  from  the  viscous  pipe 
flow  equation  and  isentropic  nozzle  equation,  one  obtains  a  value  for 
C  =  7.14  for  nitrogen.  It  would  seem  reasonable  to  assume  that  in  a 
rapidly  converging  channel  there  is  more  dissipation  and  hence  the  value 
of  C  is  greater.  The  data  inFig.Al  indicates  that  11.4  may  be  a  good 
value  for  C.  This  is  2.86  times  B  for  the  same  nozzle.  A  check  of  two 
sets  of  data  in  the  literature  on  vs.  (for  incompressible  flow, 

however)  indicates  that  this  ratio  may  be  quite  universal  (see  Fig.A2). 

If  this  is  the  case,  then  one  may  describe  the  discharge  characteristics 
of  a  given  nozzle  in  terms  of  a  single  constant  which  depends  on  geometry: 


Cd  = 


1  + 


for  Cd  ^  0.65;  Cd 


y^ed 


y^ed 

2.86B 


for  Cd  0.65  (A8) 


Since  the  Reynolds  number  is  based  properly  on  the  actual  mass 
flow,  it  is  necessary  to  do  some  rearranging  in  order  to  find  the  actual 
mass  flow  in  terms  of  that  given  by  inviscid  theory.  This  procedure  is 
shown  below. 
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^  =  6.44B2H.2d  ^d  ^  0.65  .  (A12) 


The  inviscid  mass  flow  above  which  one  uses  Eq.  (AlO)  is  from  Eq.  (A12) 


xhi  =  4.18  b2  |j,2  d 

^d  =  .65 


(A13) 
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Substitution  of  suitable  expressions  for  [I2  and  in  terms  of  y,  A2. 
A2/A1,  P]^,  P2,  and  Tg  into  Eqs.  (AlO),  (A12),  and  (Al 3)  then  yields 
general  expressions  for  compressible,  viscous  flow  through  a  nozzle 
(or  venturi  meter)  in  terms  of  quantities  which  may  be  measured  directly. 
The  results  are  given  below: 
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C.  MODIFICATION  FOR  INCREASED  KNUDSEN  NUMBER 


In  the  preceding  section,  the  effect  of  reduced  Reynolds  number 
on  the  mass  flow  through  a  nozzle  with  no  upstream  separation  was 
considered.  It  was  indicated  that  at  very  low  Reynolds  numbers  the 
discharge  coefficient  varies  as  No  lower  limit  was  given  for 

C(j>  implying  that  C^  — ►  0  as  Rg  \|^  __  0.  For  a  particiilar  liquid, 
this  probably  is  true  since  one  can  reduce  Re^j  only  by  reducing  u  or  d; 
neither  action  affects  the  basic  assumptions.  For  a  gas,  however,  one 
may  also  reduce  the  Reynolds  number  by  reducing  the  density.  This 
tends  to  increase  molecular  mean  free  path  lengths  in  the  gas.  At  a 
sufficiently  low  density  the  mean  free  path  lengths  is  a  significant 
function  of  the  nozzle  diameter  and  the  flow  may  no  longer  be  regarded 
as  a  continuum.  Deviation  from  continuum  theory  results  have  been 
found  to  occur  when  this  dimensionless  ratio,  known  as  the  Knudsen 
Number,  exceeds  0.01.  For  Knudsen  numbers  larger  than  3.0,  the  flow 
can  be  analyzed  by  free  molecule  theory.  At  intermediate  values  of 
Knudsen  number,  however,  there  is  as  yet  no  satisfactory  theoretical 
description  of  the  flow  characteristics. 

Whereas  the  continuum  theory  indicates  that  - -  0  as 

Rgjj  — ►  0,  free  molecule  theory  indicates  that  -  a  constant  0 

as  Rg^  - -  0  if  Red  is  reduced  by  reducing  the  gas  density.  This 

constant  is  a  function  of  the  nozzle  geometry  and  has  not  been  calculated 
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except  in  the  case  of  an  orifice  for  which  it  has  the  value  0.697  (Ref.  15). 
One  may,  however,  estimate  the  value  for  a  nozzle  from  the  theoretical 
results  for  free  molecule  flow  in  a  pipe  (Ref.  25).  Talking  the  pipe  length 
equal  to  its  diameter  and  P2/P1  =  0  gives 


Cd. 


min 


0.28 


The  actual  value  may  be  more  or  less  depending  upon  the  entry  geometry 
and  pressure  ratio.  Gradual  entry  shapes  would  be  expected  to  have  lower 


values  of 


For  the  nozzle  shown  in  Fig.Al,  probably 


about  0.10  when  Pz/^l  =  0.528.  One  would  expect  that  the  data  would 
begin  to  show  evidence  of  the  rarefaction  effect  for  P^  <  lOOjjiHg,  the 
pressures  for  which  the  Knudsen  number  is  greater  than  0.01.  It  is  to  be 
expected  that  will  be  reached  when  Pj  <  Ip  Hg.  For  nozzles 

with  smaller  throats,  these  rarefaction  effects  will  be  encountered  at 
higher  pressures.  For  example,  if  the  throat  diameter  is  0.01"  (0.25  mm), 
one  would  expect  that  deviations  from  continuum  results  would  begin  for 
Pj  <  4  mm  Hg. 

Although  there  is  at  present  no  experimental  verification,  it  seems 
reasonable  to  expect  that  one  could  account  for  rarefaction  effects  by 
adding  to  Eq.  (A7)  the  term 


/FX\ 

E  tanh  (  j  / 


(A17) 


where  E  and  F  are  constants  depending  on  the  nozzle  geometry.  For  the 
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nozzle  shown  in  Fig.Al,  values  of  0.1  and  3,  respectively,  seem  to  be 
appropriate.  However,  since  Eq.  (A17)  is  largely  speculative,  it  does 
not  appear  worthwhile  to  include  this  effect  in  the  derivation  of  Eq.  (A15). 
It  would  be  well,  nevertheless,  to  add  the  restriction  to  Eq.  (A15)  that  it 
be  used  only  for  X/d  <  0.01. 


D.  MODIFICATIONS  FOR  CHANGES  IN  INLET  GEOMETRY 


As  indicated  above,  the  one -dimensional  flow  equations  are  based 
on  the  assumption  that  pressure  gradients  across  the  flow  at  the  station  of 
measurement  are  negligible.  This  is  approximately  true  for  gradual 
entry  nozzles.  For  abrupt-entry  nozzles  with  short  throats  this  is 
patently  false.  It  is  common  practice,  however,  to  continue  the  use  of 
the  one -dimensional  equations  and  to  account  for  this  effect  by  means  of 
an  expansion  factor,  a  part  of  the  gross  discharge  coefficient.  The 
expansion  factor  is  that  part  which  depends  on  pressure  ratio  but  not  on 
Reynolds  number.  It  can  be  written 


Y  = 


Cd 


(A18) 


viscous 


Because  of  the  difficulty  in  obtaining  exact  theoretical  solutions 
for  flow  in  channels  which  permit  changes  in  C^j  with  pressure  ratio,  Y 
is  generally  determined  empirically.  A  theoretical  solution  was  recently 


obtained  for  inviscid  flow  through  a  two-dimensional  slit,  and  the  results 


were  found  to  agree  closely  with  experimental  data  on  flow  through  a 


sharp-edged  orifice  at  high  Reynolds  numbers.  The  interesting  feature 
of  the  theoretical  solution  (Liepmann,  Ref.  17,  discusses  this  at  some 
length)  is  that  the  maximum  mass  flow  is  not  reached  until  Pz/l’l  ^  0.05. 
This  result  is  a  consequence  of  the  fact  that  in  the  orifice,  the  centerline 
pressure  is  greater  than  that  near  the  edge.  The  centerline  flow  has, 
therefore,  a  higher  signal  speed.  Not  until  the  edge  flow  has  accelerated 
through  a  45°  angle  does  the  signal  from  it  propagate  rapidly  enough 
toward  the  centerline  to  halt  the  influence  of  lowered  downstream  pressure 
on  the  flow.  Perry  (Ref.  18)  found  that  he  could  fit  his  experimental 
results  with  the  relations 


(A19) 


for  Pg/Pi  ^  0.528 


and 


Y 


[ 


0.848  +  0.455 


Pi 


for  P2/P1  ^  0-528 


(A20) 


if  d/D  <  0.2  and  ^dy^g^-Q^g  “  i-0.  One  would  expect  some  changes  in 
Y  if  d/D  >  0.2,  but  no  data  exists  on  this  question. 
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Examination  of  Liepmann's  results  (Ref.  17)  indicates  that  for 
P2/P1  =  0.001  and  d/D  «  1,  the  variation  in  Cd..y.j^gj,Q^g  with  Reynolds 
number  is  about  the  same  as  given  in  the  handbooks  for  P2/P1  RJ  1.0. 

It  would  seem,  therefore,  that  the  viscous  and  expansion  effects  can  be 
treated  separately.  Liepmann  plans  additional  tests  to  cast  more  light 
on  this  point,  however. 

For  the  present,  if  one  assumes  that  such  separation  is  possible, 
one  may  make  use  of  the  vast  quantity  of  data  available  in  the  literature 
on  the  discharge  characteristics  of  various  geometries  measured  using 
water  or  air  with  Ap/Pj  <<  1.  From  these  data,  it  may  be  concluded 
that  a  relation  of  the  type  given  by  Eq.  (A5)  will  be  valid  for  R/d  >  0.25. 
When  R/d  <  0.25,  apparently  the  flow  will  separate  from  the  wall  at 

high  Reynolds  numbers  and  - *  0.61  as  Re^j  - -  •  Ref.  19 

indicates  that  the  expression 

Cdn  =  0.61  +  1.9  R/d  (A21) 

maybe  used  for  R/d  <  0.1.  When  d/D  >  0.2,  C^jp  increases. 

•^e  =  00 

For  example,  when  d/D  =  .4,  =  .62;  when  d/D  =  .6, 

tt0  -  00 

=  0.65;  when  d/D  =  0.8,  =  0.78.  The  dependence 

=  00  ^e  =  00 

of  C(i  on  diameter  ratio  persists  down  to  a  Reynolds  number  of  about  10, 
below  which  the  curves  are  identical.  These  results  apply  when  the 
throat  is  shorter  than  2  throat  diameters.  For  longer  throats, the  viscous 


effects  and  the  expansion  effects  are  somewhat  different,  but  they  will 
not  be  discussed  further  because  such  geometries  have  less  application 
to  the  present  study. 

It  is,  of  course,  possible  to  fit  empirical  formulae  to  each  of  these 
discharge  curves  or  to  the  family  as  a  whole.  For  example,  when 
d/D  <  0.2  and  R/d  =  0, 


Cd  =' 


I^ed  (^Gd 


1950  Red  +  2.75  (Re^  -  100)2 


(A22) 


fits  the  experimental  data  genercilly  to  within  4%  or  less.  This  relation, 
however,  does  not  lend  itself  to  simplification  or  rearrangement  as  did 
Eq.  {A5).  Even  if  one  were  to  attempt  a  two  segment  fit  such  as  was  done 
with  Eq.  (A8),  inclusion  of  dicimeter  ratio  and  leading  edge  ratio  effects 
would  hopelessly  complicate  the  expressions.  For  this  reason,  it  would 
appear  that  for  simplicity  one  should  rely  on  graphical  presentations  of 
the  characteristics  of  a  particular  device  to  determine  the  mass  flow 
represented  by  a  given  P^,  P2>  a-ud  Tg  when  R/d  <  0.25. 


E.  MODIFICATIONS  FOR  EXCESSIVE  UNFAVORABLE  PRESSURE 
GRADIENTS  DOWNSTREAM  OF  THE  THROAT 

When  a  nozzle  diverges  rapidly  downstream  of  the  throat,  a  large 
unfavorable  pressure  gradient  must  be  overcome  if  a  subsonic  flow  is  to 
adhere  to  the  walls.  If  the  gradient  is  too  large,  the  flow  will  separate 
from  the  nozzle  wall.  The  separation  point  depends  upon  the  particular 
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geometry.  Instantaneous  increases  in  nozzle  area  fix  the  separation 
point  for  all  Reynolds  numbers  and  thus  do  not  modify  the  results  obtained 
in  Subsection  A.  For  certain  relations  of  the  expansion  angle  to  the 
throat  diameter  and  length,  however,  the  wall  exerts  an  influence  on 
the  separation  point.  As  a  result,  the  separation  point  and  the  general 
flow  pattern  can  be  expected  to  respond  to  changes  in  Reynolds  number. 
The  functional  dependence  of  m  on  Re  is  therefore  quite  different  from 
that  treated  in  Subsection  A.  There  are,  of  course,  many  possible 
variations  in  geometry  which  could  be  considered;  only  one  will  be 
considered  here.  Because  of  its  relation  to  the  flow  through  screens  - 
and  possibly,  therefore,  also  to  flow  through  cloth  -  the  flow  studied  will 
be  that  over  circular  cylinders. 

From  the  drag  characteristics  of  a  cylinder,  it  appears  that  as 
the  Reynolds  number  is  increased  from  very  low  values,  the  separation 
point  moves  further  from  the  stagnation  point  to  about  the  70°  location, 
reaching  this  at  a  Reynolds  number  of  about  10^.  It  remains  there  for 
all  higher  Reynolds  numbers  for  which  the  flow  is  laminar.  This  shift 
in  separation  point  is  accompanied  by  a  decrease  in  the  width  of  the  wake. 

If  one  were  to  assume  a  nozzle  to  be  made  from  the  region  between 
two  noninteracting  cylinders,  one  can  determine  the  discharge  character¬ 
istics  from  the  experimental  drag  coefficients  of  a  single  cylinder  as 
indicated  in  Section  II.  To  perform  a  quantitative  calculation,  three 
additional  assumptions  are  required: 
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(1)  The  downstream  pressure  in  the  region  where  the  flow  is 
again  uniform  is  the  same  as  or  at  least  related  in  some  explicit  manner 
to  the  pressure  in  the  minimum  area. 

(2)  The  pressure, velocity, and  density  in  the  minimum  area  are 
those  which  would  exist  if  the  flow  were  one -dimensional. 

(3)  The  reference  Mach  number  and  Reynolds  number  are  based 
on  these  flow  conditions. 

The  force  on  an  evenly-spaced  group  of  cylinders  lying  in  a  plane 
normal  to  the  direction  of  flow  in  a  constant  area  channel  is  then 


But 


and 


AAP  =  Cq  nq2  d'.^ 


n  d'.^  =  sA 


(A23) 


{A24) 
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If  one  takes  the  drag  coefficients  for  M  ^  1.0,  Rg  >  10^  given 
in  the  literature,  one  may  compute  m/A  for  various  AP's  and  P2'®- 
the  results  are  graphed,  it  will  be  seen  that  for  Rg  ^  10^  the  curve  has 

the  same  appearance  as  that  for  an  isentropic  nozzle.  In  view  of  the 
close  similarity  between  Eqs.(A27)  and  {A14),  this  is  not  surprising; 
particularly  since  variation  in  Cq  with  M  appears  to  have  the  proper 
form  to  account  for  the  difference.  This  leads  one  to  speculate  that  the 
flow  in  region  between  two  cylinders  can  be  represented  by  the  isentropic 
nozzle  function  multiplied  by  a  discharge  coefficient  derived  from  the 
inverse  square  root  of  the  drag  coefficient  of  a  cylinder: 


as  a  function  of  Reynolds  number  is  shown  in  Fig.  A3.  For  compar¬ 
ison  purposes,  the  for  a  unseparated  nozzle  and  an  orifice  are  also 
included.  It  will  be  noted  that  the  change  in  slope  of  Cjj  vs.  Rg  extends 
over  several  orders  of  magnitude  more  for  the  two  cylinder  nozzle  than 
it  does  for  an  unseparated  nozzle  or  an  orifice.  This  situation,  reflecting 
as  it  does  changes  in  a  separated  flow,  makes  it  very  difficult,  if  not 
impossible,  to  devise  a  semitheoretical  description  in  the  manner  of 
Eqs.  (AS).  One  must,  therefore,  use  the  curve  in  Fig.  A3  in  computing 
m.  It  is  possible,  however,  to  approximate  this  curve  with  the  functions. 
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Cd  =  0.31  Re°’^^ 

Re 

< 

10 

(A29) 

Cd  =  0.43  Re°-^22 

10 

< 

Re  <  10^ 

(A30) 

Cd  =  1.0 

Re 

> 

10^ 

(A31) 

with  fair  accuracy.  Since 

Re  = . -  (A32) 

(1-s)A^2 

in  this  case,  the  actual  mass  flow  can  be  written  in  terms  of  the  ideal 
as  follows 


m 


0.31 rhi  /  d' 


([1-s]  Ay.^) 


.33 


1.5 


Re  <  10 


(A33) 


m  = 


0.43  m,  ^ \ 

U1  -s3  Ay.^) 


0.122 


1.22 


10  <  Re  <  103 


m  =  mi 


Re>  10^ 


(A34) 

(A35) 


It  is  to  be  expected  that  these  relations  will  be  valid  only  in  the 

limit  as  s  — -  0.  As  s  increases,  interaction  should  also  increase, 

particularly  at  very  low  Reynolds  numbers.  It  is  interesting  to  compare, 

1  /3 

however,  the  coefficient  inEq.  (5)  which  is  6/Re  '  with 
( l/(.43)Re' =  5.4/Re'^^^  from  Eq.  (A30).  This  indicates  that  for 
Rg  >  60,  the  interaction  effects  are  quite  small.  On  the  other  hand,  the 

right  hand  side  of  Eq.  (6)  is  considerably  different  from  (lO/Rg^/^)  ®/(l-s)^ 
indicating  that  for  Re  <  20  interaction  effects  make  the  use  of  the  two- 

cylinder  model  less  useful. 


74 


F.  CONCLUDING  REMARKS 

It  has  been  shown  above  that  it  is  possible  to  include,  without 
undue  complication  and  with  good  accuracy,  the  effects  of  viscosity 
in  the  expressions  for  flow  through  a  smooth  entry  nozzle.  For  those 
geometries  which  produce  flow  separation  (i.e.,  sharp-edged  orifices) 
no  general  analytical  description  is  known,  nor  are  the  empirical 
equations  sufficiently  simple  for  common  use.  The  effects  of  reduced 
gas  density  are  at  present  also  poorly  understood  and  little  experimental 
data  exists  to  ameliorate  the  condition.  Thus,  the  conclusion  may  be 
drawn  that  while  it  is  possible  to  account  for  some  of  the  deviations 
from  the  ideal  in  nozzle  flow,  reliance  must  continue  for  the  most  part 
on  calibration  under  actual  test  conditions. 
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FIGURE  A3 
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r  CURVES  ARE  CALCULATED  BY  EQUATION  A  26  USING 
THE  Cd  VALUES  OF  FIGURE  A3. 

2.  AREA,  SOLIDITY,  AND  WIRE  DIAMETER  OF  SCREEN 
ARE  TAKEN  AS  THE  SAME  AS  THOSE  FOR  TESTED 
SAMPLE,  15.5  cm^,  52%  AND  3.5  X  lO^^cm 
RESPECTIVELY. 
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FIGURE  5.  LOW  PRESSURE  HIGH  TEMPERATURE  PERMEABILITY  /■ 
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TEST  SAMPLE 
INSTRUMENT  PANEL 

VENTURI  METER 
INSTRUMENT  PANEL 


NOMENCLATURE 

1 

ASHCROFT  0-30  PS  1  TEST  GAUGE 

z 

54-101  •F  THERMOMETER 

3 

^HNCH  MANOSTAT  FLOWRATER 

4 

(/9-iNCH  FISCHER  AND  PORTER  FLOWRATER 

5 

VENT  COCK 

6 

OIL  MANOMETER  DISCONNECT 

r 

MANOMETER  U-TUBE 

8 

O-IOOMM  WALLACE  AND  TIERNAN  GAUGE 

9 

SILICONE  OIL  U-TUBE 

lO 

GAUGE  SEAL  OFF 

1 1 

MCLEOD  GAUGE  SEAL  OFF  VALVE 

12 

0-800  MM  WALLACE  AND  TIERNAN  GAUGE  (1  MM  DIVJ 

13 

MCLEOD  GAUGE 

t4 

i!*-inch  fine  metering  valve 

15 

2-INCH  COARSE  METERING  VALVE 

16 

WATER  OUTLET 

17 

COOLING  WATER  INLET 

18 

VENTURI  METER 

>9 

ACCESS  FLANGE 

PRESSURE  HIGH  TEMPERATURE  PERMEABILITY  APPARATUS 

I 


FIGURE  7  APPARATUS  DETAII^ 
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PERMEABILITY  OF  EXPERIMENTAL  STAINLESS 
STEEL  FABRIC 

TEMPERATUREiaSOO^K  SAMPLE  AREA  =  15.5  cm* 
DATA  UNCORRECTED  FOR  APPROACH  VELOCITY 
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PERMEABILITY  OF  MIL-C  -802IB  TYPE  II  FABRIC 

SAMPLE.  AREA^^IS.SO  cm‘ 

TEMPERATURE ~  300° K 

DATA  CORRECTED  FOR  APPROACH  VELOCITY  NON-UNIFORMITY 
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FIGURE  31 
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COMPARISON  OF  RESULTS  OF 
SEVERAL  INVESTIGATIONS  OF 
PERMEABILITY  OF 
MIL-C-802IB  TYPE  n  FABRIC 
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FIGURE  32  MIL-C-802IB  TYPE II  FABRIC  AFTER  TESTS  DURING 
WHICH  AIRFLOW  TEMPERATURE  EXCEEDED  MELTING 
TEMPERATURE  OF  MATERIAL,  DIAMETER  OF  MELTED 
SECTION  IS  ABOUT  I"  . 
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FIGURE  33  HT-I  FABRIC  AFTER  TESTS  DURING  WHICH  AIRFLOW 
TEMPERATURE  EXCEEDED  MELTING  TEMPERATURE 
OF  MATERIAL. 
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